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SECTION  I 


INTRODUCTION 


Free  flight  aerodynamic  characteristics  have  been  determined  by  several  means  of  data 
acquisition.  These  include  the  ballistic  spark  range,  wind  tunnel,  and  position  radar  coupled  with 
a  Sonde  sun  sensor  (measures  angular  orientation  relative  to  the  sun,  Reference  1 ).  The /ballistic 
spark  range  has  proven  to  be  an  accurate  means  of  obtaining  aerodynamic  characteristics. 
Reference  2,  3,  4.  The  dynamic  data  consists  of  spark  photographs  of  the  model  in  flight  at  a 
specified  number  of  surveyed  stations.  Each  station  photographs  the  model  in  both  the  vertical 
and  horizontal  planes,  thereby  allowing  for  measurement  of  the  six-degree-of-freedom  orientation 
of  the  model.  Dynamic  data  experimental  resolution  obtainable  from  a  ballistic  spark  range  is  of 
the  order  of: 

2  X  10'^  seconds  in  time 

0.1  degree  in  pitch  and  sideslip  angle 

1 .0  degree  in  roll  angle 

0.001  foot  in  X,  y,  z  coordinates 

CORRELATION  OF  FREE  FLIGHT  DATA  -  BACKGROUND 

The  problems  associated  with  direct  correlation  of  the  model  (equations  of  motion)  to  the 
test  data  are  threefold,  viz.  the  basic  differential  equations  of  motion  are  nonlinear,  the 
aerodynamic  effects  are  also  nonlinear,  and  the  test  data  are  acquired  in  a  coordinate  system 
which  is  generally  not  consistent  with  the  model  derivation. 

The  most  prevalent  method  of  analyzing  ballistic  spark  range  data  is  based  on  the  linear 
approximation  known  as  linear  theory.  Stated  briefly,  the  method  uses  a  linearized  model;  i.e.,  a 
closed  form  approximate  solution  to  an  exact  set  of  differential  equations  of  motion.  A  least 
squares  fitting  technique  is  used  to  determine  the  aerodynamic  parameters  used  in  the  linearized 
model.  Murphy  (Reference  2,  3,  4),  MacAIlister  (Reference  5  and  6),  and  others  at  (BRL) 
Ballistic  Research  Laboratory  as  well  as  Nicolaides  (Reference  6-9)  and  Eikenberry  (Reference 
10)  have  developed  this  method  to  the  extent  that  certain  types  of  nonlinearities  can  also  be 
analyzed.  However,  unless  many  cycles  of  data  are  present  and/or  multiple  experiments 
conducted,  there  is  only  one  nonlinear  aerodynamic  parameter  which  can  be  eonsisteniK 
identified,  i.e.,  the  static  overturning  moment. 
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In  analyzing  free  flight  data,  in  1969  Chapman  and  Kirk  of  NASA  Ames,  California 
(Reference  1  1)  documented  a  technique  which  allows  the  nonlinear  differential  equations  of 
motion  to  be  used  directly  in  the  data  correlation  process.  The  technique  eliminates  the 
requirement  for  closed  form  approximations  to  the  eejuations  of  motion.  It  is  essentially  a 
differential  corrections  process  of  forming  a  system  of  partial  derivatives  based  on  a  truncated 
Taylor’s  series  approximation.  The  equations  of  motion  and  partial  derivatives  are  numerically 
integrated  and  solved  for  the  aerodynamic  coefficients,  based  on  a  least  squares  formulation,  to 
best  fit  the  data.  Independent  contributions  to  this  numerical  extraction  technique  were  made  by 
Knadler  (Reference  1  2)  and  Goodman  (Reference  13),  although  the  method  is  usually  referred  to 
as  the  Chapman-Kirk  technique.  Since  this  technique  includes  a  very  generalized  model  concept, 
it  is  necessary  for  successful  application  to  have  more  insight  as  to  which  aerodynamic 
characteristics  are  pertinent  in  a  given  case.  The  previous  limitations  caused  by  nonlinearities  in 
the  equations  of  motion  and  aerodynamics  can  be  circumvented  using  this  technique. 

In  late  19()9  Whyte  and  Beliveau  (Reference  14-16)  documented  implementation  of  the 
Chapman-Kirk  technique  for  the  reduction  of  experimental  ballistic  spark  range  data.  Nonlinear 
aerodynamics  were  determined  from  a  single  data  set.  Since  1970,  the  contractor  has  been  very 
successful  in  data  analysis  of  a  variety  of  configurations  from  spark  range  data,  radar  data,  and 
Sonde  data  utilizing  this  technique  (Reference  17-22).  A  measure  of  success  might  be  the  degrees 
to  which  the  motion  data  can  be  reproduced  given  the  determined  aerodynamics  and  model.  The 
criteria  in  free  fiight  aerodynamic  data  analysis  should  be  that  the  reduction  is  not  successful 
until  the  motion  data  can  be  reproduced  to  a  probable  error  equivalent  to  that  of  the 
experimental  data. 

Until  recently,  data  reduction  investigations  on  free  flight  aerodynamic  data  have  involved 
primarily  symmetric  configurations  with  the  exception  of  wind  tunnel  investigations.  Interest  has 
been  shown  in  using  the  ballistic  spark  range  for  testing  nonsymmetric  shapes  in  free  Right  (e.g., 
elliptic  cross  section  re-entry  vehicles).  Development  of  a  data  analysis  model  and  correlation 
techniq  le  capable  of  handling  mass  and  configuration  asymmetries  was  initiated  in  January  1974 
by  the  contractor  under  government  contract.  Elimination  of  the  assumption  of  symmetry  was 
complicated  by  an  additional  coordinate  system  transformation  and  the  complexity  of  the 
aerodynamics. 

SCOPE  OF  REPORT 

The  scope  can  be  summarized  as  follows: 

•  Model  Derivation 
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Data  Correlation  Technique 


•  Test  Case  Correlation 

The  derivation  of  the  model  (equations  of  motion)  is  done  in  a  coordinate  system  rigidly 
attached  to  the  body  and  rolling  with  it.  With  the  assumption  of  symmetry,  the  model  could  be 
defined  in  a  non-rolling  coordinate  system.  Since  the  dynamic  data  are  acquired  in  an  earth-fixed 
system,  accurate  transformations  are  required  to  correlate  the  data  with  computed  values  which 
makes  the  degrees  of  freedom  highly  coupled.  Six  degrees  of  freedom  are  utilized  in  defining  the 
acceleration  and  angular  momentum  equations.  The  derivation  assumes  a  nonsymmetric 
configuration.  The  forces  and  moments  acting  on  this  configuration  are  defined  using 
aerodynamic  coefficients  to  model  the  flow  phenomenon.  Detailed  descriptions  of  the  fluid 
dynamic  effects  are  not  discussed.  Nonlinearities  in  the  aerodynamics  are  modeled  as  polynomial 
functions  of  the  components  of  the  angle  of  attack.  A  thorough  understanding  of  the  mode! 
derivation  is  imperative  to  the  overall  data  analysis  technique  due  to  the  complexity  of  the 
coordinate  system  transformations  and  modeling  of  the  flow  phenomena.  Data  analysis  using  a 
generalized  model  containing  a  large  number  of  aerodynamic  coefficients,  becomes  very  abstract 
without  this  understanding.  Physical  significance  must  be  interpreted  from  the  determined 
aerodynamic  coefficients  and  can  only  be  achieved  from  the  model  derivation. 

Numerical  integration  is  utilized  in  correlating  the  dynamic  data  to  the  model.  Both  the 
equations  of  motion  and  parametric  equations  are  numerically  integrated  using  a  fourth  order 
Runge-Kutta  Method.  As  with  the  case  of  data  analysis  of  symmetric  configurations,  the 
reduction  is  done  in  two  phases  using  least-squares  theory.  The  translational  motion  is  reduced 
first  to  determine  the  force  coefficients  and  velocity  components  required  in  the  reduction  of  the 
angular  motion.  Phase  II  embodies  the  angular  motion  reduction  to  determine  the  aerodynamic 
moment  coefficients.  Decoupling  of  the  translation  and  angular  motion  is  required  in 
implementing  the  Chapman-Kirk  technique.  Differential  corrections  using  least-squares  theory  is 
not  amenable  to  simultaneous  correlation  of  the  two  motions.  This  two-phase  data  analysis 
concept  has  been  very  successful  in  reducing  a  large  number  of  data  sets  on  symmetric 
configurations.  Reference  14-22. 

Preliminary  verification  of  the  model  derivation  and  data  correlation  technique  was  done 
using  an  ideal  test  case.  Free  flight  spark  range  data  on  the  motion  of  two  elliptic  bodies  with 
non-equal  transverse  inertias  were  then  analyzed.  The  aerodynamic  forces  and  moments  on  the 
elliptic  bodies  were  determined  with  the  resulting  probable  error  of  fit  approximately  equivalent 
to  that  of  the  measured  data. 
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SECTION  II 


MODEL  DERIVATION 

COORDINATE  SYSTEM  DEFINITIONS 

The  equations  of  motion  will  be  derived  in  a  body  fixed  coordinate  system  where  the  x-axis 
is  aligned  with  the  longitudinal  axis  of  the  missile,  and  the  coordinate  system  is  free  to  rotate 
about  that  axis  at  some  spin  rate,  In  other  words,  the  coordinate  system  is  rigidly  attached  to 
the  projectile.  The  inertial  frame  of  reference  is  the  earth.  It  is  assumed  that  the  earth  is  fixed  in 
space,  nonrotating,  and  flat. 

As  an  additional  frame  of  reference,  a  fixed  plane  coordinate  system  will  be  defined.  This 
system  is  similar  to  the  body-fixed  coordinate  system  with  the  exception  that  it  docs  not  rotate 
about  the  missile  x-axis,  and  the  y-axis  is  constrained  to  be  parallel  to  the  XY-plane  of  the  earth 
fixed  coordinate  system. 


iy,  jg,  kg  -  Earth  Coordinate  System 
'mf’ Jmf’  *^mf  “  Fixed  Plane  Coordinate  System 
'm’jm’  *^m  ~  Body-Fixed  Coordinate  System 

To  relate  the  orientation  of  one  coordinate  system  relative  to  another,  the  use  of  Euler 
angles  (0pp^  i//pn)  will  be  employed  as  illustrated  in  Figure  I,  which  relates  the  orientation  of 
fixed  plane  system  to  the  earth  coordinate  system.  The  angles  are  obtained  by  rotating  a 
coordinate  system,  initially  coincident  with  the  earth  system,  about  the  Z-axis  through  an  angle 
i//pp  and  then  about  the  Y'-axis  through  the  angle  0pp.  This  results  in  a  y-axis  which  lies  in  the 
XY-plane  which  is  part  of  the  definition  of  the  fixed  plane  coordinate  system.  To  extend  this  to 
the  body-fixed  system,  a  rotation  about  the  x-axis  is  required  as  illustrated  in  Figure  2.  The 
definition  of  the  missile  velocity  vector  components  is  also  indicated. 

EQUATIONS  OF  MOTION 

Tlie  aerodynamic  forces  and  moments  are  to  be  defined  in  the  body-fixed  coordinate 
system.  Therefore,  tlie  equations  of  motion  must  be  derived  based  on  the  acceleration  and 
angular  momentum  defined  in  the  body-fixed  system  relati\c  to  the  inertial  frame  of  reference. 
Beginning  with  Newton’s  second  Law,  the  acceleration  and  angular  momentum  will  be  defined  in 
the  body-fixed  coordinate  system,  where  tlie  center  of  gravity  of  the  piojectile  need  not  be  on 
the  x-axis. 
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X,  Y.  Z  —  Harth-Fixecl  Coordinates 

X  Y  Z  '  -  Intermediate  Coordinate  System  (after  first  rotation  -  \//) 
X,  y,  z  -  Fixed  Plane  Coordinate  System 


Figure  I .  F'ixcd  Plane  Coordinate  System 


(1) 


cl„,  V  r7  , 

F  =  111  [  +  tj  Y  X  V  ] 

L  =  ~  >•  CO  q'  X  J 

dt  ' 


(2) 


where: 


V  =  i'RB'm  +  v,^im  +  '^RB''m 

(3) 

^  T  =  P^n  +  %n  +  ^"^m 

(4) 

II 

H 

(5) 

[1]  is  the  mass  moment  of  inertia  tensor. 

The  moment  of  inertia  tensor  will  be  defined  such  that  the  body  fixed  x'y'-  plane  is  a  plane 
of  mirror  symmetry  (i.e.  “  *yz  “ 


Equations  (1)  through  (5)  will  be  stated  in  terms  of  the  body-fixed  components.  The 
following  equations  of  motion  are  obtained  after  the  vector  multiplications  and  cross  products. 


[11  = 


*x  ^xy 

~*xy  V 

-0  0 


0 

0 

U 


(6) 


Fx  =  ni(uRB  +  qwRg  -  rvRg) 

(7a) 

Fy  =  m(vRB  +  ruj^g  -  pwj^g) 

(7b) 

Fy  =  m(wRg  +  pvgg  -  qugg) 

(7c) 

Lp  =  ‘x  P  -  *xy  q  +  h  ‘1*’  +  *xy  ''P  "  *y 

(8a) 

Lq  =  ly  q  -  Ixy  P  +  lx  !’>■  -  'xy  -  'z  P*^ 

(8b) 

^r "  *  “  'xy  P^  ’y  P‘1  "  'x  P‘l  ^  'xy  ^1“ 

(8c) 

The  equations  of  motion  are  desired  in  the  form  of  linear  and  angular  acceleration  (body 
fixed)  as  a  function  of  the  aerodynamic  forces  and  moments.  Simultaneous  solution  of  Equations 
(8a)  and  (8b)  and  rearrangement  of  the  remaining  Equations  (7)  and  (8)  equations  will  result  in 
the  desired  form. 


vrb  =  PWRB  '■^'RB  + 


(9b) 


WRB  =  9URB-pvRB  +  F^/m 

,  *y^p *xy*“q  ~  *y  “  *xy  E'' +  (Ixy*"^  *y  ^*y  l/llhr 

p  = - - - 

^*x*y  “  *xy  ^ 

.  •xLq  +  lxyLp  +  '‘x  +  *y-'z)>xv9‘-  +  (lx('z-*x)  •xy“^P'‘ 

q  = - - - 

^*x*y  “  *xy  ) 


(9c) 

(10a) 

(10b) 


Lr  +  Ixy<E“  -q“)  +  (lx  -  'y)  PM 


(10c) 


Equations  (9)  and  (10)  represent  the  six-degree-of-freedom  differential  equations  of  motion 
derived  in  the  body-fixed  coordinate  system.  Once  the  definition  of  the  forces  and  moments  is 
made,  the  solution  to  these  equations  will  define  the  six-degree-of-freedom  free  llight  motion  of  a 
nonsymmetric  missile  in  body-fixed  coordinates.  The  equations  of  motion  have  been  derived  in 
the  body-fixed  coordinate  system  for  the  purpose  of  dv  fining  the  forces  and  moments.  Their 
solution  results  in  the  translational  and  angular  motion  of  the  missile  relative  to  the  body-fixed 
coordinate  system.  An  additional  set  of  transformation  differential  equations  will  be  used  for 
defining  the  motion  relative  to  the  earth.  The  data  acquired  from  a  ballistic  spark  range  is  in 
terms  of  this  motion.  Therefore,  the  transformation  equations  are  necessary  for  correlation  of  the 
equations  of  motion  to  the  data.  These  transformation  equations  are  expressed  in  terms  of  the 
body-Eixed  velocity  components  (U[^0,  Vj^g,  p,  q,  and  r),  the  fixed  plane  Euler  angles  (Opp, 
i//pp),  and  the  angle  of  rotation  about  the  missile  axis  (0). 


—  i*RB  COS0pp  COSi// pp  +  V (SI N(?ppSIN(/>C  OSi// p  p  COS0S1N i//pp) 
+  wrb  (SlN(?ppC OS0C  OSi//pp  +  S1N0SIN i//p'p) 

V Y  ~  UbbE'OS® ppSIN i// pj)  +  vgB  (SI N0 ppSlN0SIN i//p-p  +  C  OSi//p j>C  OS0) 

+  wgB  (SINOp|jCOS0SIN)|I/pp  C  OS0ppSIN0) 

Vz~  uj^B^lNOp|>  +  VBB^lN0C OSO p;p  +  "^RB^ OS0COSO j7p 


(lla) 

(llb) 

(lie) 
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0  =  P  +  TAN  0  |.|>  1  q  SIN0  +  r  COSy  1 

0  =  q  COS0  r  SIN0 
0  =  [  q  S1N0  +  r  COS0  1  /  COS(’pp 


(lid) 

(lie) 
(110 


aerodynamic  forces  and  moments 


The  rorces  and  moments  resnUin,  Irom  Ruid  How  phenom  nn.  on  nsy  >1 

rr  •  fnrm  IS  is  done  in  elassical  aerodynamics,  R. 
be  defined  m  coeificient  forn  •  R^,,-,,re„ce  23.  Dependinji  on  tlu  toptriies 

definitions  of  the  aerodynamics  are  use  ^  „ioments  conld  be  nonlinear 

and  complexity  of  the 

functions  ot  the  tota  an«  e^  ^  ^  ,s 

a,  rt,  and  (3  are  detined  in  Figure  .  .  ^  ■  iv.t<'d  with  the  exception  of  the  axial 

„,a.i,e,y  s„,aU.  Mad,  n.„n„cr  cfr.c.s  ;:  '  I  ,o  a„a,v»  f- 

,  , .. 

...■  or .0.,,  a„dc  or  a..ac.  o,  ..  ooo,oone„,s.  ..pcncncc  1,as 
l,w„  tins  nu.,l,od  or  n.oclolins  to  l,e  a.lo„uatc  (Rclccu-c  14-..). 


BASIC  FORCES  AND  MOMENTS 


The  primary  Ta'aTptsrrdisWbuti™  around  a  rotat.ng  body, 

pressure  which  is  not  coincident  with  the  center  ol  graviiy 
primary  forces  and  moments  is  as  tollows: 

C  -  Axial  Force 

-  Acts  along  the  body  axis  opposite  to  the  11^3  velocity  vector 
C  -  No 'mal  Force  (a) 


Acts  in 


the  a-pitch  plane  perpendicular  to  the  body  axis 


C, 


Normal  Force  (P) 


Acts  in  the  /3-sideslip  plane  perpendicular  to  the  body  axis 


C,p  -  Magnus  Force  (P) 


Acts  perpendicular  to  the  |3-sideslip  plane 


q 


I 

I 


Cyp  -  Magnus  Force  (a) 

Acts  perpendicular  to  the  a-pitch  plane 

Cg  -  Roll  moment 

Acts  about  x'-axis 

-  Pitching  Moment  (a) 

Acts  about  y'-axis 

C,^  -  Pitching  Moment  (/?) 

Acts  about  -  z'-axis 
Cmp-  Magnus  Moment  (|3) 

Acts  about  y'-axis 

-  Magnus  Moment  (a) 

Acts  about  z'-axis 


The  magnus  and  spin  forces  and  moments  arc  dependent  on  the  nondimensional  spin  parameter. 
pd/2v,  since  tliey  are  only  associated  with  flow  phenomena  about  a  spinning  body.  The  forces 
and  moments  will  be  modeled  as  coefficient  derivative  functions  of  the  sine  of  the  components  of 
the  total  angle  of  attack,  oT. 

'^RB 

SINa  =  -^  (12a) 


^RB 


SIN11  = 


(12b) 


Therefore: 


'^RB  ^  ''RB 

C,  =  qA(C,„  +  Cx„2<-^)  ^  ^xM  -  MRFF)  )  (13a) 


WRB  '^RB 

<^z  =  MA(q„(— )  +  C,^3(  — )  ) 


(13b) 


V 


V 


3 


''RB  ^Rli 


ty  qA(Cy^(  ^  )  +  (y|]3(  ^  )  ) 

(13c) 

pd  vj^jj 

)(  y  II 

(13d) 

pd  w,^l3 

Cyp  q  A  (  CypQ  ^  2V  *  ^  V  ^  ^ 

(I3e) 

pd 

Cg  =  q  Ad  (  Cgp  (—  )  ) 

(130 

'^RB  3 

C,„  =  q  Ad  (CV^(— )  +  C„^^3(  ^  )  ) 

(13g) 

'"RB  vrb  3 

Cn  ""  4  Ad  (  ^ (— — )  ) 

(13h) 

pd  Vr3 

^mp  ^  4  Ad  (  C,„p^  (— )  (-^ )  ) 

(13i) 

pd  Wr,j 

C|,p-4  Ad(C„p^(— )(  ^  )) 

(13j) 

DAMPING  MOMENT 


Tlie  yaw  clamping  moment  being  considered  is  defined  as  acting  perpendicular  to  the  body 
axis  but  independent  of  the  yaw  (angle  of  attack)  plane  orientation.  It  is  dependent  on  the 
angular  velocity  components,  q  and  r.  Modeling  of  the  damping  moment  is  accomplished  using 
Uie  aerodynamic  coefficients  and  The  sense  of  these  components  is  depicted  in  Figure 

qd  w,^3  “ 

Cmq  =  BAcl(— )  ) 


^  nr  B  Ad  (  ^^  )  (  )  > 


(1 4b) 


It  should  be  noted  that  a  more  exact  definition  of  the  damping  moment  would  be  to 
consider  it  as  being  dependent  on  the  cro.ss-angular  velocity  components  q'  and  r'  which  are 
coincident  with,  and  normal  to.  the  angle  of  attack  plane,  respectively  (Reference  23).  These 

components  are  shown  in  Figure  4.  The  more  exact  definition  would  introduce  crosstlamping 
terms  into  the  equations  of  motion. 
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BODY-FIXED  ASYMMETRIES 


Body-fixed  aerodynamic  forces  and  moments  are  due  to  misalignment,  cant,  or  asymmetry 
of  body  and/or  lifting  surfaces.  Figure  4  depicts  a  center  of  gravity  which  is  completely 
nonplanar  relative  to  the  primary  body  axes  (x'-y'-z').  These  forces  and  moments  are  commonly 
referred  to  as  trims.  The  existence  of  trims  would  cause  the  free  fiight  angular  motion  of  a  body 
to  oscillate  about  an  axis  which  is  not  coincident  with  the  initial  body  x'-axis.  The  body  fixed 
forces  and  moments  will  be  defined  by  the  aerodynamic  coefficients  Cy^,  C^y,  and  Cy^y  as 
illustrated  in  Figure  4.  The  forces  Cyy  and  C^y  act  at  the  center  of  pressure  causing  the  moments 
Cyy  and  C,yy  ubout  the  center  of  gravity,  respectively.  Trims  along  and  about  the  body  X'-axis 
are  absorbed  in  the  axial  force  and  roll  moment  coefficients  as  defined  in  liquations  (13a  and  f). 


AERODYNAMIC  EFFECT  DUE  TO  CROSS  FLOW  ORIENTATION 


These  effects  are  due  to  the  orientation  of  the  aerodynamic  surfaces  with  respect  to  the 
cross  flow.  The  orientation  of  these  surfaces  will  be  defined  by  the  aerodynamic  roll  angle,  0', 
which  is  the  clockwise  rotation  by  the  aerodynamic  surfaces  with  respect  to  the  cross  flow  as 
illustrated  in  Figure  5.  The  cross  flow  velocity  vector  is  the  resultant  velocity  perpendicular  to 
the  primary  body  axis,  x'. 


0'  =  TAN‘ 


''RB 
1  _ 

wrb 


-r 


(15) 


The  phase  angle,  f,  due  to  the  initial  orientation  of  the  plane  of  symmetry  relative  to  the 
body-fixed  coordinate  system  axes  should  typically  be  zero.  The  most  significant  induced 
aerodynamic  effects  resulting  from  missile  orientation  relative  to  cross  flow  are  generally  the 
induced  rolling  moment,  the  induced  side  force,  and  the  induced  side  moment  (Reference  23). 
Since  the  cross  flow  will  have  some  effect  on  normal  force  and  pitching  moment,  they  will  be 
considered  also.  These  forces  and  moments  are  a  function  of  the  total  angle  of  attack  in  addition 
to  the  aerodynamic  roll  angle. 


SIN  2  5  SIN  ( .10'  ) 

=  Induced  side  force  acts 

(1 6a) 

perpendicular  to  yaw  plane 

'n0^SIN2  5siN(n0') 

=  Induced  side  moment 

(1 6b) 

'N0^S1n2  5S1N(,i0') 

=  Induced  normal  foice  acts  in 
yaw  plane 

(1 6c) 

WaS'N"aSIN(n0') 

=  Induced  pitching  moment 

(I6d) 
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CROSS  FLOW 
VELOCITY  VECTOR 


HORIZONTAL 


PLANE  OF  SYMMETRY 


Figure  5.  Aerodynamic  Surfaces  Relative  to  Cross  Flow 


=  Induced  rolling  moment 


(I6e) 


Cq0^,SIN  -  a  SIN  (  110'  ) 

n  =  number  of  axially  symmetric  fins  or  surfaces 

GRAVITATIONAL  FORCE 

The  force  of  gravity  will  be  cielmed  in  vector  notation  as  follows: 

Gravity  Force  =  -  mg  ( 1 7) 

Replacing  the  unit  vector  by  its  representation  in  the  body  fixed  coordinate  system,  the 
gravity  force  will  be  resolved  into  components  which  are  consistent  with  the  definition  of  the 
aerodynamic  forces. 

-mgkj.  =  mgSIN0ppi„^  -  mgCOSdppSIN0j,„  -  mgCOSOppCOS0'k„,  (18) 

It  is  assumed  that  the  earth  is  flat,  and  the  earth’s  rotation  may  be  neglected. 

SUMMARY  OF  FORCES  AND  MOMENTS 

Combining  the  forces  and  moments  which  have  been  defined  in  Kquations  (12)  through 
(18),  the  terms  Fy,  F^,  Lp,  L^j,  will  be  defined  for  substitution  into  Fquations  (9)  and 
(10). 

'^RB  "  vrP  2 

Fx  =  qA[Cxo  +  ^'xa2<— 7"^  +  ^'x^2  (M  -  MRFF)  ]  +  mgSINOpp  (19) 

''RB  ''RB 

Fy  =  OA  I  So  ^  "  S’/33  (— )  +  Spa 

'"RB  (20) 

-CN0aS«N2  a  SIN(n0')(— )  +  Cy^^SIN^  a  SIN(n0')(— ^)  ] 

-  mgCOS0p’pSlN0 

wr3  w,^3  2  pjl 

S  =  M  A  [  CVo  +  (— )  +  (— )  +  C,p^  ^ 

-V  V 

,  '^RB  ''Rli 

“^N0a‘''"^“  ®  SIN(n0')(— ^)  Cy^^SIN^  a  SIN(n0')(— ^)  ] 

mgCO,SOpp('OS0 
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pd 

Lp  =  qAd(  Ccp(-:^)+  Cg0^SIN-aSIN(n0') 


'^RB  '^RB  ^ 

-  q  Ad  [  C  (  ^)  +  C  (  —  )  +  (":^K  ^ 

pd  Vrb  Wr,^ 

+  (— )  +  <^tn0aS*N2  a  SIN(n<&')(— ^) 

''RB 

+  <^'n0a^'N2aSIN(n0')(— )] 


'^RB  ''RB  ^  ''RB  ~ 

Lr  =  nAd(C„„  +  C„^(— )  +  C„^3(— )  +  +  C„r2  (— >  > 


pd  wru  vr3 

'*'  ^npa  "  ^m.AaSIN2  a  SIN(n0')(— ) 


+  C„^S1N- a  SIN(n0')(— ^)  J 


SI  CTION  III 


DATA  CORRI  LATION  TliCHNIQUi; 

DIFFERENTIAL  CORRECTIONS  (SENSITIVITY  EQUATIONS) 

The  essential  steps  of  the  analysis  a.e  briefly  summarized.  The  exact  equations  of  motion 
are  derived  containing  all  necessary  aerodynamic  coefficients.  By  integration  of  the  equations  of 
motion  the  theoretical  trajectory  is  obtained.  The  actual  trajectory  is  known  from  the  ballistic 
range  test  data.  The  gist  of  the  method  is  to  adjust  the  aerodynamic  coefficients  contained  in  the 
equations  of  motion  in  such  a  way  that  the  generated  trajectory  matches  the  actual.  The 
mathematical  treatment  is  a  sensitivity  analysis.  The  sensitivity  equations  or  parametric 
differential  equations  are  ordinary  differential  equations  whose  dependent  variables  are  the 
partial  derivatives  of  the  state  variables  with  respect  to  the  parameters.  The  sensitivity  equations 
are  integrated  in  parallel  with  the  equations  of  motion  to  yield  sensitivity  coefficients  (partial 
derivatives),  which  reflect  the  sensitivity  of  the  computed  solution  with  respect  to  each 
aerodynamic  coefficient  or  system  parameter  (e.g.,  initial  condition). 

The  use  of  a  truncated  Taylo.’s  series  expansion  brings  about  a  method  of  qiiasilineariz.ation 
which  uses  these  sensitivity  coefficients  to  linearize  the  change  in  the  solution  of  the  nonlinear 
equations  of  motion  due  to  a  change  in  the  aerodynamic  coefficients  or  system  parameters.  As  an 
illustration  we  will  consider  the  following  functional  representation  of  the  solution  to  a  system  of 
dynamic  equations. 

Let: 


n  ■*  o  ** 

0^  =  tf^,  t)  =  computed  solution  to  the  0  equation  of  motion  based  on  the  initial  estimates 
to  the  system  parameters,~c°. 

c  =  (ci^  ^’2'  *'3 . ‘■pi*  “the  system  parameter  vector  containing  parameters  to  be 

determined. 


^c 


=  g(c°,  t)  =  computed  solution  to  the  i//  equation  of  motion  based  on  initial  estimates. 
=  hfir^,  t)  =  computed  solution  to  the  0  equation  of  motion  based  on  initial  estimates. 
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0|;(t),  0|;(t)  =  experimental  data  at  discrete  points  in  time. 


Hmploying  a  Taylor’s  series  expansion  about  what  will  be  defined  as  tlie  nominal  solution,  and 
neglecting  second  and  higher  order  terms  results  in  the  following: 


_  _  P 

0,,  =  f  (C°+  AC,  t)  =  f  {C'^,  t)  +  S  ( - )  A  C: 

i=i  ac|  ' 

(25) 

_  _  P 

=  g  (CO+ AC,  t)  =  g  (CO  t)  +  2( - )  AC: 

^  i=i  ac]  ' 

(26) 

_  p  90^0 

0.,  =  h  (C0+  AC,  t)  =  h  (CO,  t)+  S  ( - )  AC: 

i=l  9C| 

(27) 

Where  P  is  the  number  of  system  parameters  to  be  determined,  liquations  (25)  through  (27) 

represent  the  nominal  solution  which  provides  the  optimum  correlation  to  the 
Data  acquisition  errors  or  measurement  noise  is  defined  in  liquations  (28). 

experimental  data. 

Rq  =  0p  (t)  -0^ 

(28a) 

R^  =  i^[.;(t)-0^. 

(28b) 

R0  =  0b  (t)  -  0^. 

(28c) 

The  assumption  made  is  that  in  the  absence  of  measurement  noise,  the  dynamic  equations  of 
motion  exactly  define  the  experimental  data  based  on  the  system  parameters,T.  Inclusion  of  the 
measurement  noise  in  liquations  (25  through  27)  gives  the  following: 


p  dO^o 

Rfl  =  0 11  -  f  (C®  +  AC,  t )  =  0 f  (C®,  t )  S  ( - )  AC: 

i=l  dC\  ' 

^  P  3i//j,o 

Rx  =  I//,;  g(C«+ AC,  t)=  -  g(C^  t)-  V(__)  AC 

i=l  I 


(29a) 


(29b) 


p  d0„o 


R^  =  0n  -  h  (C^^  +  AC.  t)  =  <t>y.  -  h  (C^,  t)  -  - )  AC: 

0  L  t  gj..  I 


(29c) 
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Let: 


i;o  =  o,;-i  (c«  t) 
L^  =  .^n-g(c«  t) 
L^  =  0,..  -  I.  (CO,  t) 


(30a) 


(30b) 


(30c) 


Theoretically  R^,  and  R^  have  been  defined  as  the  measurement  noise.  In  an  iterative  sense, 
R  is  one  iteration  ahead  of  E,  which  is  defined  in  Equations  (30).  In  the  limit  as  convergence  is 
achieved  E  should  approach  R  being  a  measure  of  the  measurement  noise.  E'^,  Ej^,  and  E^  can  be 
thought  of  as  an  approximation  of  the  measurement  noise  based  on  initial  estimatesreO. 

LEAST-SQUARES  THEORY 


A  suitable  technique  for  matching  the  theoretical  trajectory  to  the  actual  is  least-squares 
theory.  Through  successive  iterations,  the  residual  functions,  R^,  R^,  and  R^,are  minimized  to 
the  measurement  accuracy  of  the  experimental  data.  The  residuals  to  be  minimized  are  the 
differences  between  the  experimental  data  and  the  computed  solution  which  are  defined  in  the 
previous  section  as  the  measurement  noise  E.  Application  of  the  theory  results  in  the 
minimization  of  the  squares  of  the  residuals  based  on  variations  (Ac)  of  the  system  parameters 
(c).  The  square  of  the  residual  functions  to  be  minimized  are  as  follows: 

N  _  _  N  p  dO^o  ^ 

R0“  =  Z  [  Op-  -  f  (CO  -t-  AC,  t)  1“  =  Z  1  Ep  -  Z( - )  ACj  )  (31a) 

k=l  ^  k=l  i=l  ^ 


N 

R^2  =  i^(  ^  g(CO+Xc,t)]2 
k=l  '  k 


N  p 
k=I  ^  i=l 


1: 


N 

R  “  =  Z  [(Ap  -  h  (C° 
^  k=l 


N 


p  d<p„o 


+  AC,t)|-  =  Z(E.-Z(-— )  ACj)- 


^  k=l 


i=l 


(31b) 


(31c) 


Minimization  of  the  sum  of  the  squares  of  the  residuals  will  be  achieved  by  taking  the  partial 
derivative  with  re.spect  to  Ac  and  setting  it  equal  to  zero.  This  is  indicated  in  Equations  (32).  N  is 
defined  as  the  total  number  of  data  points. 
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Equation  (33)  may  be  solved  to  provide  eorrections  ([AC])  to  be  applied  to  tlie 
aerodynamic  coefficients  and  system  parameters  (c^),  knowing  tlie  partial  derivatives 
([A]  —  sensitivity  coefficients)  and  the  approximated  measurement  noise  ([E]-e.g. 

APPLICATION  TO  EQUATIONS  OF  MOTION 

In  correlating  ballistic  spark  range  data  using  least-squares  theory,  analysis  of  the  swerve 
motion  (x,  y,  z)  is  decoupled  from  the  angular  motion  {0,  i//,  <j>).  The  equations  of  motion  (EOM) 
used  in  the  decoupled  swerve  analysis  are  as  follows: 

SWERVE  EOM 

PAW- 

)[  -  CxCQSOppCOSi//p|,  Cy  (SINdppCOS.//ppSIN0  SINi//ppCOS0)  (34a) 


q  (SIN0ppCOSi//ppCOS0  +  SIN;//ppSIN(A) 

_  Pd 

+  Cyp  (~)(SlNOppSIN0COS^pp  SlNi//ppCOS0) 
_  pd 

-  qp  (— )(SlNOppCOS0COS^pp  +  SIN0SIN0PP  )  ) 


PAW- 

Vy  =  (-^^)  I  -  C^COSOppSIN^pp  -Cy  {  SINOppS!N0ppSlN0  +  COS0ppCOS0  )  (34b) 


-  Cy  (SINOppSIN0ppC  OS0  ^  COS0ppSIN0) 

_  pd 

+  Cyp  (— )(SINOppSIN0ppSIN0  +  COS0ppCOS0) 

_  pd 

Cyp  (— )(SINOppSIN0ppCOS0  -  COS0ppSlN0)  ] 
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(34c) 


PAV- 

V7=(-^^)|  Fusing, FySiN^cosop,, 

C'^  COS0C  OSO|;|,  +  Cyp  (— )  SIN0COSO,r|, 

_  Pd 

-  C,p  {— )  COS0COSO,;|,  I  g 

Equations  (34)  were  derived  by  differentiating  the  transformation  Equations  (11a),  (11b),  and 
(11c),  and  substituting  Equations  (9a),  (9b),  and  (9c)  for  U|^g,  V|^g,  and  Wj^g.  Definition  of  the 
aerodynamics  in  Equations  (34)  follows. 

_  wgg  2  vj^g  2 

=  ^xo  +  <^xa2  >  +<^x/J2<— )  +C\^m-MR\IF)  (35a) 

'^RB  '^RB  ^ 

Cy  =  Cy„  +  Cy^  (— )  +  Cy^3  (  — )  (35b) 

_  '^RB  '^RB  ^ 

=  C/.o  +  <^’za  ^za3  7">  (35c) 


'^RB 

C  =r  ( - 

yp  ^ypa^  Y 

_  ''RB 

('x.p "  ^  zp(3 


(35d) 

(35e) 


Direct  correlation  of  experimental  data  to  the  equations  of  motion  for  the  determination  of 
the  aerodynamic  coefficients  requires  a  system  of  partial  differential  equations,  sometimes  called 
sensitivity  equations.  This  system  of  equations  results  from  taking  partial  derivatives  with  respect 
to  each  aerodynamic  coefficient  on  the  equations  of  the  motion. 
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The  system  of  partials  required  for  direct  correlation  to  the  swerve  or  position  data  (x,  y,  z) 
is  as  follows: 


[  A]'  = 


ax 

ax 

ax 

ac, 

dC2 

aCr 

av 

av 

av 

aci 

ac2 

aCf 

az 

az 

az 

aC| 

ac2 

aCr 

(36) 


These  partial  derivatives  are  obtained  from  numerical  integration  ofdWJdC-  d\ JdC:,  and 

avacj. 

For  simplicity  of  derivation  as  well  as  computational  purposes,  partial  derivatives  will  be 
stated  in  generalized  form.  Derivation  of  the  partial  derivative  for  Equation  (34a)  follows. 


3v„  av  aCi 

—  =  A,—  +K,j— ^ 
aCj  acj  •'  acj 


(37) 


where: 

PAV  w,^l3  -  v,^g  - 

Cxa2(— )  +  2C  x^2(— )  +  (3M  -  2^! REF)  )  COSO ppCOS,//,: 


''RB  ^RB 

~  2(CyQ  +  Cy|3  (— ^  )  +  Cy^2  )  )  ( S 1 N 0 p |)COS p p S l N ^  —  SlNi//pprOS0) 
'^RB  ^RB  ^ 

~  ^^^zo  ^  ^za  ^  ^  ^  ^za3  ^  y  ^  ^  (SINOppCOS0ppCOS0  +  SIN0ppSIN0) 

wrb  pd 

Spa  ^  (SIN0ppSIN0(  OS0pp  SIN0ppCOS0) 

vrb  pd 

~  S.p|3  ^  Y  ^  (  2y  ^  (SIN0ppCOS0COS0pp  +  SIN0SlNi//pp  )  J 
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3 


pAV“ 

c.p.  K|= - (SINOp,,COS0ppCOS0  +  SlN\i/p,,SIN0) 

for  j  corresponding  to 

'>2  2  1/'’ 

V  =  (V,  +Vy  -HV,  ) 

dv  av^  av^  v,,  av., 

_  ^  + _ y  y  +  *  _ _ 

aCj  V  aCj  V  acj  v  acj 

As  indicated  in  the  preceding  partial  derivative,  those  variables  computed  during  analysis  of  the 
angular  motion  are  assumed  known  and  therefore,  are  not  treated  as  dependent  variables.  The 
variables  which  are  assumed  known  are  W|^g^Y>  ''RB/V*  ^FP’  '^FP'  Derivation  of  the 

complete  set  of  partial  derivatives  used  for  analysis  of  the  position  data  is  contained  in  Appendix 
A. 


Th;  equations  of  motion  used  in  the  angular  motion  analysis  are  as  follows: 

ANGULAR  EOM 

lyLp  +  l^yLq  -  (Ijj  +  ly  -  \y)  IxyP'^  +  dxy"  V  ^•y  " 

(38a) 

,  *xki  *xy*"p  ^*x  ^  *y  -  ^j)  *xy  “  *x^  ~  *xy"^ 

q  - - ^ - 

^*x*y  “  *xy"^ 

(38b) 

Lr  + Ixy  (p^  -  q^)  +  (lx  -  lylpq 
r  - 

(38c) 

vrB  ""  P'^RB  ~  '■“RB  ^  SlN0COSOp|, 

(38d) 

^RB  ""  ‘P'RB  “  P'^RB  ^  ^  COS0('OSO|.  |, 

t38e) 

0  =  p  +  TANOj.p  ((]  S1N0  +  r  COS0) 

(38f) 

0  =  q  COS0  r  SIN0 

(38g) 

•  • 

The  0  equation  is  not  included  since  0p|>and  0  are  not  needed  lor  the  analysis. 
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Definition  of  the  aerodynamics  in  Equations  (38)  follows. 


pd 

Lp  =  qAdl  Ces  +  Cyp  SlN2aSlN(n0')l 


(39a) 


L„  =  qAd  I  ■*  <^’nia  ^  v  ^  ^  ^ 


3  qil  qd 


W{^[j  '''RB 


pd  v,^u  '^RB 

+  C-mpS  ^  “  SIN  (no—  I 


''RB 

+<^n0a  SlN2aSIN(n0')(— )1 


VRB  ''RB 

Lr  =  qAd  [  ^’n|3  ^ni33  ^  v~ 


rd 

■•■  ^  nr  ^  2V~^ 


rd  V|^3 


pd  W(^i^ 


+  Cnr2  ^  ^  2V  V 


''RB 

)  C,„^aS'N-aSlN(n0')(— ) 


(39c) 


"'RB 

+  C„^aS‘N-aSlN(n0')(— )1 


''RB 


''RB 


pd  W|^ij 

)  ■*■  ^'vnn  ^  -ix;  \/  ^ 


F Y  =  q A  I  Cy o  +  Cy^  (— )  +  Cyp3  (— '  '-y pa  s V  V 


''RB 


w 


RB 


CN0aSlN“  a  SlN(n0')  (-^)  +  Cy^a^^N*"  « 


(39d) 


Fz  “  qA  1^7.0  ^  ^  za  ^  Y 


WRB  '^RB  ^ 

)  +  ^za3  ^  2V  V  ^ 


(39e) 


"'RB 


''RB 


_CN0a^lN2aSlN(n0')(— ^)-<^'y0a^'N‘'aSlN(n0')(  ^  )1 
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The  system  of  partial  differentia!  equations  required  for  direct  correlation  to  the  angular 
motion  data  "^RB’  follows; 


^^RB 

3^RB 

^^RB 

9C‘i 

9(’-, 

3«//rb 

3\^rb 

dii'RB 

9('l 

9C2 

«’s 

. 

90 

9Cl 

DCS 

These  partials  are  obtained  from  numerical  integration  of  9v|^g/9Cj,  9wj^g/9Cj,  90/9Cj. 
Having  obtained  3vj^g/9Cj  and  9w  rb/SCj,  the  transformation  equations  are  used  to  obtain 
partials  of  and  \I/gg  as  follows: 


*RI)  =  -  SIN-' (  vkb/V) 

Orb  =  SIN-1  (wrb/(VCOSi)-rb)  ) 
Using  the  chain  differentiation  rule  results  in  the  following: 


^®RB 

^'’RB 

90RB  Dwrb 

^^RB 

9Cj 

3C3 

9vrb  3Cj 

9wrb 

3'//rB 

^''RB 

90RB  3wrb 

1 

^'^'RB 

9Cj 

-9C3 

9vrb  3Cj 

9wrb 

The  derivatives  of  Equations  (41 )  are; 

^  =  0.0 
9wrb 

9i//rb  _  1 

avRB  (v2-vrb^) 

^^RB  _  _ • _ 

9wrb  (v2cOs2\(/rb- wrb“^ 
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(41a) 

(4Ib) 


(42a) 


(42b) 


(43a) 


(43b) 


(43c) 


% 


I 


ovrb  (V^COS^^I^B  -  wrb^) 


wrb  S1Ni//rb 

COS1//RB  -  vrb“) 


(43cl) 


Equations  (43)  are  substituted  into  Equations  (42)  to  obtain  ^Orb/^^J  ‘'”^1  •"'e 

used  to  form  the  matrix  [B],  Equation  (40).  The  system  of  correction  equations  used  to  correlate 
the  angular  motion  data  is  stated  as  follows  in  matrix  notation. 

(Bl^  [W1[B]  1  AC]  =  [  BjT  1  W)1  R  ] 


1  0  0 
(W]=  0  1  0 

0  0  W3 


Equation  (44)  represents  a  weighted  least  squares  formulation  to  correlate  Orb-  '/^rB' 

[W]  matrix  is  used  to  weight  the  correlation  of  0  since  the  measurement  errors  on  0  are  of  the 
order  of  5  to  10  times  larger  than  Orb  and  <|'rb- 

where: 

W3  =  weighting  factor  on  0 

Derivation  of  the  partial  derivatives  for  Equations  (38)  is  contained  in  Appendix  A.  The  velocity 
is  assumed  known  from  analysis  of  the  swerve  motion  and  is  therefore  not  treated  as  a  dependent 
variable. 

The  yaw  data  correlated  during  analysis  of  the  angular  motion.  Orb  and  ^rr.  is 
transformed  after  the  swerve  analysis  as  follows: 

''RB  =  Vx(  SIN0ppSIN0COS0p,,  -  COS0S1N0PP  ]  ,45^. 


Vyl  SlN0ppSIN0SIN0pp  +  COS0COS0PP 
V2SlN0COSOpp 


Wrb  =  Vxl  SINOppCOS0COS0pp  +  SlN0SlN0pp 


(4.Sb) 


+  VyI  SlNOppCOS0SlN0pp-  S1N0COS0PP 


+  V7COS0COSOPP 


(45c) 


'/'RB - ^  (''RB^^^ 


Orb  =  sin-« 


'^RB 

^  VCOS^RB  ^ 


(45d) 


In  summation,  the  measured  six-degree-of-freedom  data  (X,  Y,  Z,  0pp,  >//pp,  0)  are 
correlated  as  follows.  The  position  data  are  fit  to  Equations  (34)  resulting  in  a  determination  of 
the  aerodynamic  force  coefficients.  The  yaw  data  (0pp  and  \//pp)  are  then  transformed  to  the 
body-fixed  coordinate  system  using  the  computed  earth-fixed  velocity  romponents,  Vy,  and 
V^.  Equations  (45)  are  the  transformation  equations  used.  The  angular  motion  data  (Orr.  ’/'rb* 
0)  are  then  analyzed  using  Equations  (38)  to  result  in  a  determination  of  the  aerodynamic 
moment  coefficients.  A  weiglited  least  squares  formulation  is  used  to  account  for  the  larger 
measurement  errors  in  0  relative  to 
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SnCTlON  IV 


PROGRAM  UTILIZATION 


Data  acquired  from  the  Lglin  Aeroballistic  Range  Facility  consist  of  pliotographic  records 
of  the  angular  and  linear  position  of  a  body  in  free  flight.  The  pliotographs  are  taken  from 
surveyed  stations  spaced  on  the  order  of  5  to  20-foot  intervals  along  the  length  of  the  range.  The 
records  are  obtained  from  dual-axis  spark  shadowgraphs.  Optical  axes  at  each  shadowgraph 
station  are  mutually  orthogonal  with  the  range  centerline.  A  film  reader  is  used  to  result  in 
numerical  interpretations  of  the  shadowgrams  which  provides  histories  of  attitude  and  position  as 
functions  of  time  (X,  Y,  Z,  Opp,  i//pp,  <t>  vs.  time). 

The  six-degree-of-freedom  data  is  then  correlated  to  the  equations  of  motion  as  derived  in 
the  body-fixed  coordinate  system.  The  aerodynamic  model  used  is  derived  in  Section  II  and  the 
data  correlation  technique  is  detailed  in  Section  111. 

REDUCTION  PROCEDURE 

The  normal  reduction  sequence  to  determine  the  force  and  moment  coefficients  consists  of 
the  execution  of  five  computer  programs.  The  final  two  computer  programs  are  for  simultaneous 
reduction  of  up  to  three  rounds  of  similar  configurations.  Figure  6  represents  a  basic  block 
diagram  of  the  reduction  sequence  illustrating  the  flow  of  information  between  programs.  A 
description  of  the  function  of  each  program  follows. 

ROLLNUT  •  Preliminary  transformation  from  earth  fixed  to  missile  angles. 

•  Fits  missile  angles  using  modified  linear  theory. 

•  Computes  dense  profile  (relative  to  data  profile)  of  body-fixed  angles 
and  angular  rates  from  approximate  ..ansformations  for  use  in 
ROLL-HEliVE. 

•  Computes  initial  conditions  for  ROLL-llFEVE  and  ROLL-ANGLES. 


ROLL-llFFVF 


Provides  a  best  fit  to  the  position  data  using  earth  fixed  equations  of 
motion. 


Provides  a  preliminary  fit  to  the  roll  data  (0)  assuming  Iw 

Ay 


=0. 


T|  t  ■  -jcg. jg.  n  ;u  r.-  ,-^y 

»■  n  --  ■»  -jifcaMBlMditeiMn  I II iini tfnM«.>iWti*r>ii 


•  EARTH  FIXED  ANGLES  VS.  TIME 

•  TRANSLATION  COORDINATES  VS.  TIME 


REDUCTION  OF: 

•  MISSILE  ANGLES  USING 
LINEAR  THEORY 


REDUCTION  OF: 

•  TRANSLATION  DATA  WITH  EARTH 
FIXED  SWERVE  EOM 

•  ROLL  DATA  WITH  APPROX.  EOM  (lj^Y  =  0) 


DATA 


MANGLE 
DATA 
ACCUMULATOR 


I 


AERO-MOMENT 
COEFFICIENTS  . 


ACQUIRED 

RANGE 

DATA 


•  PHYSICAL  PROPERTIES  OF  MISSILE 

•  ATMOSPHERIC  CONDITION  OF  RANGE 


RESULTS;  (APPROX.) 

•  BODY-FIXED  ANGLES 

•  BODY-FIXED  ANGULAR  RATES 

•  INITIAL  CONDITIONS 


RESULTS; 

•  AERODYNAMIC  FORCE  &  ROLL  COEFFICIENTS 

•  TRANSFORMATION  TO  BODY-FIXED  ANGLES 

(EXPERIMENTAL  DATA) 


ROLL-HEEVE 
X,  Y,  Z,  0 


■AERO-FORCE 

COEFFICIENTS 


I 


ROLL-ANGLES 
^RB'  '^RB'^ 


/  MUL-HEEVE  \ 
DATA  J 

ACCUMULATOR  J 


DATA 


REDUCTION  OF; 

•  BODY-FIXED  ANGULAR  MOTION  DATA  TO  FULL  6DOF 
EOM 


RESULTS:  AERODYNAMIC  MOMENT  COEFFICIENTS 


1 


•  SIMULTANEOUS  REDUCTION 

MANGLE 

OF  SIMILAR  CONFIGURATIONS 

MUL-HEEVE 

®RB-  ’^RB-  ^ 

•  DETERMINES  OPTIMUM  COMMON 

X.  Y,  2,0 

AERODYNAMICS  FOR  CONFIGURATION 
SIMILARITIES 

I  MULTIPLE  REDUCTIONS 


Figure  6.  Body-Fixed  Data  Analysis  System 
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•  Transforms  the  fixed-plane  angle  data  (lApp,  i/'pp)  to  body-fixed 

angles  (0pjj,  using  the  computed  earth-fixed  velocity 

components  and  experimental  roll  data. 

•  The  total  velocity  profile  is  computed  for  use  in  ROLL-ANGLl'S. 

•  Results  in  the  aerodynamic  force  and  roll  moment  coefficients. 

ROLL-ANGLES  •  Provides  a  best  fit  to  the  angular  motion  data  (Orr-  '/'RB’ 

the  body-fixed  equations  of  motion. 

•  Computes  a  dense  profile  of  the  body-fixed  angles  and  angular  rates 
for  a  final  swerve  analysis  (ROLL-HEEVE). 

•  Results  in  the  aerodynamic  moment  coefficients. 

MANGLE  •  Simultaneous  angular  reduction  of  up  to  three  data  sets  from  similar 

configurations. 

•  Determines  optimum  common  aerodynamics  corresponding  to 

configuration  similarities  and  unique  coefficients  per  data  set  to 
those  parameters  which  are  not  common. 

MUL-HEEVE  •  Simultaneous  swerve  reduction  of  up  to  three  data  sets  from  similar 

configurations. 

•  Determines  optimum  common  aerodynamic  force  coefficients 

corresponding  to  configuration  similarities  and  unique  coefficients 
per  data  set  to  those  parameters  which  are  not  common. 

The  input  required  to  execute  each  program  is  defined  in  Appendix  B. 

During  implementation  of  the  derived  model  (equations  of  motion)  and  data  correlation 
technique,  options  were  included  such  that  the  data  analysis  could  be  made  using  symmetric 
aerodynamics.  These  options  exist  in  the  form  of  an  additional  set  of  aerodynamic  coefficients 
paralleling  those  derived  in  Section  11.  For  example,  the  pitching  moment  coefficients  and 
are  paralleled  by  which  is  the  symmetric  pitching  moment  term.  The  purpose  of  these 
options  is  three-fold.  If  the  data  being  analyzed  is  on  a  symmetric  configuration,  then  tlie 
symmetric  aerodynamics  should  be  computed  while  setting  the  nonsymmetric  coefficients  to 
zero.  Secondly,  if  only  slight  asymmetries  are  present  in  certain  forces  and  moments  and  the 
quality  of  the  data  is  not  good  enough  to  determine  the  asymmetric  effects,  then  the  forces  and 

moments  may  be  constrained  to  be  equal  in  each  of  the  body-fixed  planes  by  using  the 

symmetric  terms.  Finally,  when  initial  estimates  of  the  aerodynamics  are  not  very  good,  initial 
execution  of  the  programs  should  be  done  using  the  symmetric  terms  such  that  convergence  will 
be  assured  and  improved  initial  estimates  determined  for  the  final  analysis.  A  definition  of  the 
complete  set  of  aerodynamic  coefficients  and  equations  of  motion  for  each  program  follows. 


Tlic  program  ROLL-HliEVU  which  docs  the  swerve  analysis  includes  a  preliminary  analysis 
of  the  roll  data.  The  equation  of  motion  was  simplified  by  assuming  Ij^y  equal  to  zero.  This 
correlation  of  the  roll  data  is  done  independently  of  the  correlation  to  the  position  data. 

The  sign  convention  of  the  aerodynamic  coefficients  in  some  cases  differs  from  the 
derivation.  The  signs  reflected  in  the  following  equations  are  exactly  as  programmed.  The  user  is 
advised  to  use  these  equations  to  verify  the  sign  convention  of  determined  coefficients  resulting 
from  the  computer  programs. 


ROLL-NUT  I'niiations  of  Motion 


X  X 

0,„  =  K|e  ’  SIN(0| +co,X-H  W|X“/2) 

X  X 

+  K20  “  SIN  (02  +  wiX ^2x2/2) 

+  K^SIN  (0-^  +  0) 

X  X 

0,„  =  K|C  '  COS(0, +co,X  + J},X-/2) 

X  X 

+  K2C  “  cos  (02  +  CO2X  +  a32X“/2) 

+  K^^  OS  (02  +  0) 

ParaniL'tors  to  be  tletcrminocl  from  correlation: 

Kj  =  nutation  amplitude  vector,  degrees 

X|  =  nutation  damping  exponent 

0|  =  initial  orientation  of  nutation  vector,  degrees 

a)|  =  nutation  frequency,  deg/fl 

• 

CO  I  =  rate  of  clumge  of  nutation  frequency,  deg/ft- 

1^2  ~  precession  amplitude  vector,  degrees 

X2  ~  precession  damping  exponent 

02  “  initial  orientation  of  precession  vector,  degrees 

CO2  =  precession  frequency,  deg/ft 

•  ^ 

CO2  =  rate  of  change  of  precession  frequency,  deg/ft“ 

=  trim  amplitude  vector,  degrees 
03  =  initial  orientation  of  trim  vector,  degrees 
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3rr: 


ROLL-HEEVE  Equations  of  Motion 

pAV- 

~  ^ ^  ^  ~  Cj^COSflppCOSi//pp  -  Cy  (SlN0ppCOSi//ppSlN0  —  SlNi//ppCOS0) 

-  q  (SlN7?ppCOS0ppCOS0  +  SlN0ppSIN0) 

+  Cyp  (— )  (SINOppSIN0COS0pp  -  SIN0ppCOS0) 

_ 

-  ^zp  (SIN0ppCOS0COS0pp  +  SIN0SIN0pp)  ] 

pAv2  _ 

Vy  =  (— — )  [  -  CxCOSOppSlN0pp  -  Cy  (SlNOppSlN0ppSIN0  +  COS0ppCOS0) 

-  q  (SINOppSlN0ppCOS0  -  COS0ppSlN0) 

_  pel 

+  Cyp  (— )  (SlNOppSIN0ppSIN0  +  COS0ppCOS0) 

_ 

-  C^p  (— )  (SlN0ppSIN0ppCOS0  -  COS0ppSlN0) ) 


pAv2 

[  -  CASING  pp  -  CySlN0COS0pp 

_  _ 

-  C2COS0COS0pp  +  Cyp  (— )  SIN0COS0pp 
pd 

-Czp(-^)COS0COSOppl-g 

.  pd  ly  _ 

P  =  t  Cfip  +  Cg5  +  C^  SIn2  a  SIN(n0')  ]  +  ( — j - )  qr 

X  I Y 


0  =  p  +  TANOpp  (q  SIN0  +  r  COS0) 


ROLL-HEEVI:  Equations  of  Motion  -  Continued 


n  =  no.  of  fins,  i.e.,  no.  of  symmetries 


''RB 

0'  =  TAN'^  ( - ) 

wrb 


wrb  '"rb  r  9- 

Cx  =  C,o  +  q^2(— )  )  +^xM<M-MREF)+,C,-„2^1N  a 


vrb  ''RB  '^RB  ''RB  I 

Cy  =  Cyo  +  Cy^  (— )  +  Cy^3  ^  -"I^Na  CNa3SIN  a  (— )  ^ 


WRB  ^RB  ^1  '^RB  '^RB 

Cz  =  C,o  +  (— )  +  C,^3  (— )  +1Cn5(— )+  ^ 


WRB  I  '^RB  * 

Cyp  ^  ^ypa  ^~v~^  ^  V  ^  I 


Vrb  I  '"RB  ’ 
^zp  ^zpp  "^l^ypa  ^  V  ^  I 


-  SYMMETRY  TERMS 


ROl  L-ANC’il,l'S  l'x|iiations  of  Motion  (  ontinucd 


wIktc; 


RU 


''RB  ^  ' 


''RB 


'RB 


Normal  Typo  ('y  “  ^  YO  “  ^ 


pci  wj^B 


pci  wrb 


Magnus  Typo  ^Vpa  ^  V 


l^OypH'V— )] 


''RB 


liulucocl  Normal"'  ^  N0a  SIN  (n0  )  S1N“  a  (  ^  ) 


w 


RB 


Inclucocl  Siclo"  +  <^'Y05  SIN  (n0')  SIN- 5  (— ^) 


'^RB 


'^RB 


'^RB 


Normal  Typo  ^Z"“^ZO  ^  Za 


t'Za3<-v - '  ;-t'Nd<— >-<^Na3S'N‘a(- 


PJ  v. 


Magnus  Typo 


'RB 

<^Zp(3  <•— ’ 


pd  vrb 


^  Ypa  ^  2V  ^  ^  V  ^ 


w 


RB 


Incluoocl  Normal’*  ^  N0a  ^  SIN"  a  <  y  ) 


'RB 


Incluoocl  Sicio"  t’Y^a  SIN  (n0')  SIN-  5  (— ^) 


"  Not  inoluclocl  in  ourront  program. 


I  SYMMI  IRYTPRMS 


I’itchini; 


Dampini; 


Magnus 

liuluccil  Pitching 
( rrec].  effect ) 

Induced  Side  Moment 
(damping  effect) 

Pitcliing 

Datnping 

Magnus 

Induced  Pitching 
(fret|.  effect) 

Induced  Side 
(damping  effect) 

Roll  Trim 
Spin  Decay 


ROLl.-ANCU.I  S  laiualions  of  Motion  Continued 
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'RB  '^RB  '  .  '^RB 

^  m  “  ^  mo  ^  ma  ( — v — ^  ^  ma3  ^  ^  ^  ma  ^  v  ^  ^  ma3  ®  * 


w 


qd 

qd 

"'RB  ■ 

1  qd 

1 

qd 

mq  ^"Ty 

'  ^  niq2  ^ 

V  ^ 

■*^1^1117}  ^”2^^  ■*■  ^mq2 

^2V 

pd  VRB  I  pd  Vrb  I 

^  ^mp^3  ^  2V  *  ^  I 


w 


RB 


+  C,„^-SIN(n0')SIN-a(— :^) 


V 
''RB 


+  Cn^5:SIN(n0’)SIN-5(~-) 


''RB  ''RB  ''RB 

^n  "  ^no  “  ^n/3  ^  y  ^  “  ^n03  ^  ^  v"'^  ~ 


RB 


rd 


rd  V 


RB 


rd 


rd 


''-I 


+  ^nr  +  <^nr2  >  +^nq  <-3^)  +  t'mq2  SIN*^  a 


pd  wrb  I  pd  wrb  I 
^npa  ^“3^^  ^  '’’[^npa  ^”3^^  ^  y  ^  1 


''RB 


V 

"'RB 


+  C’n0aSIN(n0')SIN-a(-^) 


(’.,=  C 


SYMMI  TRY  TP  RMS 


+  ('■,.  SIN‘a  SIN  (n0') 


Induceil  Roll 


SKCTION  V 


DATA  ANALYSIS  RLSULTS 


Since  the  model  derived  in  Section  II  represents  a  larj>e  class  of  Hying  objects,  it  is  necessary 
to  extract  by  coefficient  variation  that  one  which  corresponds  to  the  actual  missile  as  observed 
by  tests,  specifically  the  free  Hight  motion  of  a  nonsymmetric  body.  The  basic  form  of  the 
body-fixed  six-degree-of-freedom  equations  of  motion  is  common  to  the  field  of  Free  Flight 
Dynamics.  The  basic  aerodynamic  forces  and  moments  assumed  to  be  acting  on  a  nonsymmetric 
body  as  discussed  in  Section  II  have  been  used  in  missile  and  aircraft  design.  Modeling  of  the 
nonlinearitics  in  the  forces  and  moments  as  well  as  induced  effects  is  not  nearly  as  common  and 
is  difficult  to  determine  from  the  literature. 

Preliminary  checking  of  the  programs  described  in  the  previous  section  was  achieved 
utilizing  an  ideal  test  case  generated  using  a  six-degree-of-freedom  trajectory  program.  The  test 
case  was  essentially  that  of  a  symmetric  cone.  The  purpose  of  using  such  a  test  case  was  that 
during  development  and  implementation  of  the  programs  a  major  concern  in  verification  was  to 
insur  that  the  transformations  were  correct  in  going  to  body-fixed  coordinates,  the  differential 
corrections  procedure  was  properly  implemented,  and  in  general  the  entire  reduction  procedure 
was  consistent.  This  approach  proved  valuable  towards  preliminary  verification  in  solving 
problems  during  implementation. 

Ballistic  spark  range  data  on  two  elliptic  cones  were  provided  by  AF.DC  for  analysis.  The 
cones  had  non-equal  transverse  inertias  (ly  which  provides  an  indication  of  the  asymmetry 
involved.  The  physical  properties  are  tabulated  in  Table  1.  Correlation  of  the  free  fiight  range 
data  to  the  six-degree-of-freedom  equations  of  motion  was  done  using  the  reduction  procedure 
detailed  in  Section  IV.  The  results  of  the  analysis  are  tabulated  in  Table  2.  I'hese  aerodynamic 
coefficients  result  in  a  probable  error  of  fit  to  the  data  which  is  about  equivalent  to  the 
acquisition  errors.  The  elliptic  cone  was  treated  as  a  two  fin  configuration  for  the  induced  effects  j 

since  a  plane  of  symmetry  is  encountered  every  180  degrees.  The  maximum  yaws  for  the  two  w 

elliptic  cones  were  about  6  and  3  degrees  for  Cones  No.  1  and  2  respectively.  Figures  7  and  8 
illustrate  the  free  Hight  angular  motion  of  each  cone  as  plotted  in  the  body-fixed  coordinate 
system.  These  plots  illustrate  the  correlation  to  the  angular  motion  data  resulting  from  single  j 

reduction  fits.  Figure  9  illustrates  the  fixed  plane  angular  motion  of  Cone  No.  1 .  ! 
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TABLE  1.  PHYSICAL  PROPERTIES 


CONE 

NO. 

Vni 

(Ft/sec) 

DIAM. 

(Inches) 

LENGTH 

(Inches) 

WEIGHT 

(Pounds) 

>x 

(Lb-in^) 

•y 

(Lb-in“) 

■z 

(Lb-in") 

*xy 

(Lb-in") 

1 

8410. 

1.484 

4.704 

0.2293 

0.04158 

0.35180 

0.33933 

0. 

2 

8395. 

1.484 

4.699 

0.2338 

0.04220 

0.35640 

0.34350 

0. 
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Simultaneous  swerve  and  angular  motion  reductions  were  performed  utilizing  the  programs 
MUL-HFEVE  and  MANGLFS.  The  result,  Table  2,  of  this  final  reduction  was  a  better 
determined  set  of  common  aerodynamics  with  a  unique  set  of  trim  coefficients  for  each  cone. 
The  spin  or  roll  rate  dependent  aerodynamic  effects  (magnus  and  spin  decay)  were  not  computed 
since  the  roll  rate  was  very  small. 

The  accuracy  of  the  determined  force  and  moment  coefficients  must  be  discussed  in  terms 
of  their  importance  towards  affecting  the  observed  motion  of  each  cone.  The  angular  motion  of 
Cone  No.  1  was  about  twice  the  magnitude  of  Cone  No.  2  and  contained  a  roll  reversal  whereas 
No.  2  had  an  increasing  roll  angle  throughout  the  flight.  Intuitively,  it  should  be  expected  that 
the  correlation  of  data  for  Cone  No.  1  would  result  in  better  determined  coefficients  for  the  2- to 
6-degree  angle  of  attack  region,  and  Cone  No.  2  for  theO-  to  2-dcgree  region.  Assuming  that  both 
cones  represent  the  same  configuration,  simultaneous  correlation  to  both  sets  of  observed  motion 
should  result  in  determined  aerodynamics  covering  the  0-  to  6-degree  angle  of  attack  region. 

The  computed  trim  moments  (C^^q,  were  of  small  magnitude,  but  well  determined 
since  they  are  dominant  terms  in  the  equations  of  motion.  They  result  in  trim  angles  of  0.214 
degree  for  Cone  No.  1  and  0.291  degree  for  Cone  No.  2.  The  pitching  moment  coefficient 
derivatives  (Cm^,  Cn^)  were  well  determined  and  differed  by  about  28  and  21  percent  for  Cones 
1  and  2  respectively.  A  roll  induced  pitching  moment  term  was  computed  for  Cone  No.  1 . 
Inclusion  of  this  term  into  the  model  being  correlated  to  the  data  resulted  in  reducing  the 
probable  error  of  fit  to  the  angular  motion  data  by  10  percent.  The  induced  pitching  moment  is 
plotted  as  a  function  of  aerodynamic  roll  angle  and  angle  of  attack  in  Figures  10  and  11 
respectively.  The  contribution  of  the  induced  effect  to  the  pitching  moment  is  of  the  order  of  3 
percent  for  Cone  No.  1 .  The  correlation  to  the  angular  motion  for  Cone  No.  2  was  not  improved 
by  including  the  induced  effect  in  the  reduction.  This  is  probably  a  result  of  the  smaller  angular 
motion  and  may,  in  part,  account  for  the  difference  in  Cn^  between  the  single  reduction  fits. 
Simultaneous  analysis  of  the  two  cones  resulted  in  a  common  pitching  moment  with  induced 
effect  included. 

The  induced  roll  moment  (Cg^)  was  very  critical  in  correlating  the  roll  data  since  the  spin 
rate  was  very  small.  The  roll  history  for  Cone  No.  1  contains  a  roll  reversal.  Figures  12  and  13 
illustrate  the  induced  roll  moment  coefficient  plotted  as  a  function  of  aerodynamic  roll  angle  and 
angle  of  attack,  respectively. 

The  damping  is  probably  the  most  difficult  parameter  to  determine.  It  is  a  dynamic  term 
and  considering  the  time  period  (0.1  second)  of  the  data,  its  effect  on  the  motion  should  be 
relatively  small.  The  damping  moment  coefficient  determined  was  consliained  to  be  ecpial  in  each 
body-fixed  plane.  Computation  of  unique  moment  components  in  each  plane  resulted  in  no 
improvement  in  the  correlation  to  the  data. 
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Simultaneous  analysis  of  the  swerve  data  resulted  in  valid  axial  force  and  normal  force 
coefficients  for  the  0-  to  6-degree  angle  of  attack  region.  Table  2  indicates  variation  in  the 
coefficients  determined  from  individual  fits  to  the  data.  However,  considering  the  distribution 
and  magnitude  of  the  angular  motion  data  (Figures  7  to  9),  the  actual  forces  computed  are 
consistent  between  the  individual  and  multiple  correlations  to  the  data.  Cone  No.  2  had  small 
angular  motion.  Therefore,  computation  of  axial  force  as  a  function  of  angle  of  attack  was  not 
meaningful.  It  is  felt  that  simultaneous  analysis  provided  a  good  distribution  of  data  to  result  in 
an  accurate  computation  of  the  force  coefficients. 


SECTION  VI 


CONCLUSIONS 


SUMMARY 

A  model  (equations  of  motion)  has  been  developed  which  is  capable  of  computing  the  six- 
degree-of-freedom  trajectory  of  a  nonsymmctric  missile.  The  derivation  of  the  dynamic  equations 
does  not  include  the  effects  of  the  Earth’s  curvature  and  rotation.  The  forces  and  moments 
considered  to  be  acting  on  the  missile  are  described  in  a  body-fixed  coordinate  system.  The 
inertial  frame  of  reference  for  the  model  is  the  Earth.  Aerodynamic  coefficients  are  used  to 
define  the  forces  and  moments.  They  are  assumed  to  be  nonlinear  functions  of  the  total  angle  of 
attack,  components  of  the  angle  of  attack,  and  the  aerodynamic  roll  angle.  Eolynomial 
expansions  are  made  as  functions  of  the  sine  of  the  aforementioned  angles.  Mach  number 
dependencies  of  the  forces  and  moments  luive  not  been  included  with  the  exception  of  the  axial 
force.  Modeling  of  the  forces  and  moments  has  been  done  to  include  a  large  class  of  Hying 
objects.  The  generalized  model  is  suitable  for  application  in  the  analysis  of  ballistic  spark  range 
data. 


Ballistic  spark  range  data  provide  the  actual  trajectory.  The  analysis  of  that  data  re(iuires 
adjusting  the  aerodynamic  coefficients  of  tlie  equations  of  motion  such  that  the  theoretical 
trajectory  matches  the  actual.  A  least  squares  data  correlation  technique  has  been  developed 
which  forms  the  basis  of  a  sensitivity  analysis  u.sed  to  adjust  the  coefficients.  A  system  of 
sensitivity  equations  of  partial  differential  equations  are  developed  based  on  a  truncated  Taylor's 
series  expansion.  These  are  derived  from  the  six-degree-of-freedom  equations  of  motion.  Using 
the  least  squares  method  of  correlation,  analysis  of  the  translational  motion  is  decoupled  from 
the  analysis  of  the  angular  motion. 

Five  computer  programs  have  been  developed  to  result  in  a  data  analysis  system  for 
correlating  ballistic  spark  range  data  or  a  nonsymmctric  configuration.  The  first  program  of  the 
series  (ROLLNUT)  correlates  the  angular  motion  data  to  a  modified  linear  theory  yaw  eejuation. 
The  primary  objective  of  this  program  is  to  compute  initial  estimates  of  the  linear  aerodynamic 
coefficients  and  initial  conditions,  and  approximations  for  the  body-fixed  angles  and  angular 
rates.  The  second  program,  ROLL-MEEVE,  is  used  to  determine  the  aerodynamic  force 
coefficients  through  correlation  of  the  translational  data  (X,  Y,  Z)  to  the  earth-fixed  equations 
of  motion.  Correlation  of  the  angular  motion  data  0)  to  the  body-fixed  e(]uations  of 


motion  is  done  next  in  tlie  sequence  using  ROLL-ANCiLES  to  determine  the  aerodynamic 
moment  coefficients.  These  three  programs  are  capable  of  computing  the  optimum  set  of 
aerodynamic  force  and  moment  coefficients  to  minimize  the  probable  error  of  fit  to  a  set  of 
ballistic  spark  range  data.  The  remaining  two  programs  are  used  for  simultaneous  correlation  of 
up  to  three  sets  of  data  on  similar  configurations  to  determine  a  common  set  of  aerodynamic 
coefficients.  The  program  MUL-HEEVE  is  used  to  correlate  the  translational  motion,  and 
MAN(il/E  correlates  the  angular  motion  data.  The  equations  of  motion  and  sensitivity  equations 
are  numerically  integrated  using  a  fourth  order  Runge-Kutta  technique. 

Simulated  motion  was  used  for  preliminary  checking  of  the  model  and  data  correlation 
technique.  Ballistic  spark  range  data  on  two  elliptic  cones  was  analyzed  using  this  data  analysis 
system.  The  resulting  probable  error  of  fit  was  approximately  the  experimental  accuracy  of  the 
data.  The  side  force,  normal  force,  pitching  moment,  yawing  moment,  and  trim  forces  and 
moments  were  determined.  A  roll  in  iiic  •  moment  dependent  on  the  aerodynamic  roll  angle  and 
square  of  the  total  angle  of  attack  was  also  determined.  Simultaneous  analysis  of  the  two  data 
sets  resulted  in  a  well  determined  .set  of  common  aerodynamics  for  the  elliptic  cones. 

RECOMMENDATIONS 

As  previously  stated,  the  model  represents  a  large  class  of  flying  objects.  Therefore,  given  an 
actual  trajectory,  there  are  many  possible  combinations  of  aerodynamic  coefficients  available  to 
consider  in  determining  those  which  represent  the  forces  and  moments  on  the  missile. 
Information  contained  in  the  system  of  sensitivity  equations  could  aid  in  determining  which 
aerodynamic  coefficients  provide  the  most  accurate  representation  in  matching  the  actual 
trajectory.  This  should  be  considered  in  future  investigations. 

The  aerodynamic  model  derived  is  felt  to  be  accurate  and  complete  for  the  analysis  of 
ballistic  spark  range  data  on  a  nonsymmetric  configuration.  It  has  been  successfully  used  in  tlie 
analysis  of  two  elliptic  cones.  However,  it  should  be  noted  that  there  are  many  possible 
formulations  of  the  aerodynamic  forces  and  moments.  As  more  experience  is  gained  in  analyzing 
asymmetric  configurations,  improvements  in  the  manner  of  modeling  the  aerodynamics  may  be 
desired.  These  improvements  will  probably  be  in  the  functional  relationships  of  the  forces  and 
moments  on  the  angle  of  attack,  angle  of  attack  components,  and  aerodynamic  roll  angle. 

Least-squares  theory  has  been  used  as  the  basis  for  the  sensitivity  analysis  which  represents 
the  data  correlation  technique.  The  swerve  motion  is  decoupled  from  the  angular  motion 
resulting  in  a  two-phase  data  analysis  concept.  This  concept  has  been  very  successful  in  analyzing 
a  large  number  of  data  sets  on  symmetric  configurations  ( Reference  14-22)  as  well  as  in  analyzing 
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the  two  elliptic  cones  on  the  body-lixed  reduction  programs.  However,  the  Maximum  Likeliliood 
Method  is  also  a  suitable  data  correlation  technique.  This  technique  can  be  used  to  correlate  all 
six  degrees  of  freedom  simultaneously  for  a  determination  of  the  aerodynamic  force  and  moment 
coefficients.  Brielly,  this  method  has  the  capability  of  introducing  statistical  information 
pertaining  to  the  measurement  accuracy  into  the  reduction  process  (Reference  24).  Ballistic  spark 
range  data  has  been  analyzed  using  the  Maximum  Likelihood  Method  (Reference  25).  It  is 
recommended  that  consideration  be  given  for  utilization  of  this  technique  in  future  investigations 
towards  eliminating  the  decoupling  of  the  swerve  and  angular  motion  analysis. 
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APPHNDIX  A 


HQUATIONSOF  MOTION 

PARTIAL  DERIVATIVES  SWERVE  ANALYSIS 

Derivation  of  the  partial  derivatives  used  in  the  swerve  analysis  follows.  Parameters 
computed  during  analysis  of  the  angular  motion  are  W[^jj/V,  vj^g/V,  Opp,  i//pp,  and  0.  These  are 
not  treated  as  dependent  variables. 

The  aerodynamics  used  in  the  swerve  analysis  are  defined  as  follows: 


•  _  f  ]  r_  -  - 


CO3^ppCO50  +  S/V  Pfp  Sttv  {^) 

+■  ^yp  (  2V  )G  inSfp  S/N0CO6(f^fp-  5/^  COS  0) 

-  C_rp  C~^)  (^thteppCa50CC^%p-f-'5lN0  <,tN  (f^Fp)  j 
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Taking  the  partial  derivative  of  Vy'. 


(A-4) 


where: 


f'C?Cx,^i>Cx.a(^^^/+2C>,p2(  W^C^(3/vA-2/)lP£Fjljco56f,5/W^ 

~  2Cy  (•sif/efp"SiN<;^fp^iN0  ^- <^ff  coz>0  ) 

-  2Cj_  (,^IKlQff.SlU  tPfp  CoS^-CoS^^-fp  ) 

^  c-  (W  (  ^iNOfp  S/fv0fp  y.  CO^<Pfp  CO^0  ) 

-  (^)(^rS'GFp5lSf<0>fpCO^0-COi<^FP^^U0)  J 


'^2j  '  t-ttapec/  /c  //e  J Co-C  /c/9 

Consider  the  V,^  HOM: 


Vj  -  C  2  ^  5lfv0pp  ~  ^1^0  CO 5  Qfp 

~  Co^  0  C05  &fp  *-  Cy^  (2^ )  S/V<^CO50p| 

-  ^  )  COS0  cos  0e„  ]  " 


Taking  the  partial  of  V^: 


3  _ 


where: 


A,  ,  /r  £Cj 

nr. 


(A-6) 


^3  ^  2(n  )[  ■*  (  V  )  V  )  ^^>,m('2>M-  2/riPcF))  "5/N'6fp 

~  2C^  (5tt^^C05Qrp)- 2C^CcoZipcosefp)  *  Cy^,  H'i)  (si  u^Cc-^Opp) 

c  rp  C  2W  )  (^Co-b(^  COStOpp  )  J 

Equations  A-1,  A-3,  and  A-5  arc  the  equations  of  motion  used  to  correlate  the  position  data. 
Equations  A-2,  A-4,  and  A-6  are  the  corresponding  partial  derivatives.  A  preliminary  analysis  of 
the  roll  data  (0)  is  done  during  the  swerve  analysis.  The  following  equations  of  motion  are  used 
which  assume  that  the  product  of  inertia,  L.„  is  zero. 

Ay 


(A-7) 


0 

0  -  P  *  Opp  ( °l  Zisf  0  F  rc.o^ 


(A-8) 


wherr : 


n  -  no.  of  fins  or  symmetries 


0'  =  aerodynamic  roll  angle 


=  TAN  ■'  ( - ) 

'^RB 
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ANGULAR  MOTION  ANALYSIS 


Derivation  of  the  partial  derivatives  used  in  the  angular  motion  analysis  follows.  The  total 
velocity  is  treated  as  a  known  based  on  the  swerve  analysis. 

The  aerodynamic  forces  and  moments  used  in  the  angular  motion  analysis  arc  as  follows: 


-  *1  AJ  [^i<i  ^  (zv  )  S/v'  ^  ] 

C.r\0oi  C  V  J  J 

.  C„,  C  v’j  (  ^-J\  c..  (^)(^‘) 


(A-11) 


(A-12) 


^  N  0^J  C  y~  J 


(A-13) 


(A- 14) 


(A-15) 


Note:  The  induced  normal  and  side  forces  indicated  by  h'quations  (39a)  and  (39e)  are  not 
included  in  the  programs. 
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For  convenience  of  derivation,  the  following  identities  will  be  stated 


■5  /  v  *  “ 

5/W  0') 

PFINV  = 

PFIUV^  ~ 

P5fS/tV  = 

P5F5/4VV  = 

PIN'D  VV  - 

PJWZ)  V/VY  = 


'V'ci  i- 


V' 


=  -5  /fv 


(  n  TAn'i^  V**)  ) 


5/ wCocO 


3 


=  iTOS 


f  — — 


..  _  co.<n0';  I 

Par>,o  [  V,*  J 


3  ^-,6 


Q  SIN^^  5/ v(^0*!) 


3 


=  ^/n'S-CPF/Nw) 


«e 


?  "'f’ug 


X^t  = 

(  Ik  ~  )  ^*Y 

rs,?  = 

(iKy  -Xy(Ty-T^)) 

ir5  3  = 

(  r,r,.  -^*y) 

ZS4  = 

(Xn  C-Tz-rx)--^*/) 

Mipilllf" 


The  dependent  variables  to  be  considered  in  taking  partial  derivatives  are  p.  q,  r,  Op.-p,  0, 
Vj^l^,  and  Wj^g.  Derivation  of  the  partials  follows. 

Consider  the  P  EOM: 

P  .  Lrl^.  tbx.  P'- 

rs3 

Taking  the  generalized  partial  of  P: 


9_P  ^ 


0 


'3JL 


^  S3  ^  ^  8^ 

^  3c 


»8  +  6s  f  i 

3Cj 


^-£i 

;>C; 


(A-17) 


where; 


6, 


^^3  ^  ^  ZS3  2S3  3  6 


5, 


J-52  ^  ^  -Txy  ^ 

r53  X-33 
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rS2  ^  -T'S/ 

-r^i3  ^  2-53  ^ 


2h 

5V.J 

rs3 

-  -tLi 

^io 

i’M- 

r-53 

-^3 

Partial  of  P  with 

respect 

to  the  j  coefficient 

Consider  the  q  l  OM : 


rs3 


(A-18) 


Taking  the  generalized  partial  of  q; 


D, 

^  D.  S3  . 

-  7)3  ^ 

yc  y  Dr 

PC- 

3Cj 

3c, 

PC; 

2£ 

PC: 


(A-19) 


where: 


D.  = 


2).  = 


<:r34) 
-  r  y- 

P-  s  y.  «>^ 

X53 

X53 

3  P  rs3 

JTi 

Pis 

•TSJ 

XS3 

9*) 

JS3 

(X54j 

r.63 

P-  = 


Tid  p-v;, 
_  X  >-  1 


Jmx  'P^p 
X-S3  PV,g 


D  =  ^ 

15  3 


»s 


K^,j  =  Partial  of  q  with  respect  to  the  j  coefficient 
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Consider  the  r  EOM: 


Consider  the  wnn  l:OM: 


^ Ft /m  -  ^cos<^  cos  e, 


(A-24) 


Taking  the  generalized  partial  of 


If  -  Mr  Mr  pi  *  If  ^  ^ 


(A-25) 


where: 


G,  ^  -■\r.R^ 


G-,  - 


_  p  '  2E* 

_L  3^* 

^  Z 


~  9  Cc^0  •^/fvQf 


C-i  ~  CC^Bf, 


Consider  the  0  hOM. 


P  ^  rf\uefpi'=\ ^ 


Taking  the  generalized  partial  of  0. 


'D<^ 


2i=  w,  ||'"»3Ci  ' 

OC] 


where; 


H,  =  1-0 


/Vj  =  r  AN  Q  Fp 


<)/^0 


/y  =  rAN  Ofp  cct></' 


/y^  -  "c>FP  C=f  “S/.v  J  ^ 

/7y  -  rAtNCfr  (^  ^  ^  ^ 


Consider  the  0  EOM. 


a 


5  -=  <^cos0-rs/.vj2^ 


Taking  the  generalized  part\al  ol 


i'O: 


-jQ  ^  91  ~  ^/r^0  ^ 

9Q  =  COSJZ)  pc,.  ^  ^  3Cj 


Taking  partial  derivatives  of  liquations  (A-1  1 )  through  (A-1  5)  results  in  the  lollowing. 


^  ^  P  -  °l  Ad  2'v  J  ^  ^  p 

^  P 

?ii>  =  ^Ad  (P6FCAY) 
p/._o  =  (P^FSAWj 


Ad  (^)  (  ^) 

P 


PZc)  _ 


^  Ad  (svj  C'y  )  ) 


3A7e 

=  ‘7/iJ  [c.w  (vj  ^  "■ 

51j:  ^  9  Ac/  (2^ )  V* ) 

pp 

?_ir=  7/^c/ (if  ^ 

r  A  '  \  4-  C  ( )  ^■  Cf^rx\2'^  )\  V*/ 

2-k.  ^fij[c„f,(-;;)<-‘-‘'^^'^  V*  ^ 

?ir. 

-I 


5-^-i:  =  cj  /I 
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(Tlie  reverse  of  tliis  page  is  blank) 


RQLLNUT  INPUT 


HARD  NO.  1  (10A3) 


IBM  Col 

Variable 

Descriptor 

1-30 

SHOTN 

Title  Card 

CARD  NO. 

2  (717) 

1-7 

LN 

Start  point  of  fit 

8-14 

NA 

End  point  of  fit 

15-21 

ND 

Incremental  section  length 

22-28 

NAUTO 

0  input  W] ,  ^2 

1  input  CMAE 

29-35 

NROLL 

0  roll  data  on  card  No.  3 

2  roll  data  on  card  No.  5 

36-42 

NCOORD 

0  AF.DC  data 

1  Eglin  or  test  case  data 

43-49 

NBODY 

0  Output  in  fixed  jilane 

1  Output  in  body  axis 

Normal 

Settings  -  Eglin 

AEDC 

NAUTO 

0 

0 

NROLL 

0 

2 

NCOORD 

1 

0 

NBODY 

1 

1 

CARD  NO.  3a. z  (7F10.0) 

1-10 

TG 

Time,  seconds 

11-20 

XC 

X  ,  inclu's 

range 

21-30 

YG 

Y  ,  inches 

range 

31-40 

ZG 

Z  ,  foot 

range 
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CARD  NO.  3a 

,z  (7F10 

.  0)  (continued) 

Al-50 

DMG 

m 

direction 

51-60 

DNG 

n  cosines 

61-70 

DPE 

p  or  0  (see  NROLL) 

r  1 

CARD  NO.  A 

(IFIO.O) 

1-10 

TG 

TG  =  200  (stop  card) 

Note:  Card 

No's  5 

and  6  are  not  supplied  when  NROLL  =  0. 

CARD  NO.  5 

(2E1A.6. 

IFIA.O) 

1-lA 

PHI 

^FP 

15-28 

DUMMY 

Z  -  not  used 

S 

29-A2 

DUMMl 

CARD  NO.  6 

(2E1A.6. 

IFIA.O) 

1-lA 

PHI 

15-28 

DUMMY 

29-A2 

DUMMl 

DUMMl  =  200  (stop  card) 

CARD  NO.  7 

(7F10.0) 

1-10 

DIA 

diameter,  inches 

11-20 

AIX 

I  ,  lb  in^ 

X 

21-30 

AIY 

I  ,  lb  in^ 

y 

31-AO 

AIZ 

1  ;  lb  in 
z 

Al-50 

AIXY 

1 

xy 

51-60 

WGT 

weiglit,  lbs. 

61-70 

ACC 

CG  oil  set,  inches 
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CARD  NO.  8  (AFIO.O) 


1-10 

XCGR 

Reference  CG,  percent  of  length 

11-20 

XLR 

Reference  length,  inches 

21-30 

XCG 

CG,  percent  of  length 

31-40 

XFIN 

No.  of  assymetries  or  fins 

CARD  NO. 

9  (3F10.0) 

1-10 

RHO 

3 

density,  slugs/ft 

11-20 

ASOUND 

speed  of  sound,  ft/scc 

21-30 

VMACH 

reference  niach  no. 

CARD  NO. 

10  (4F10.0) 

1-10 

CNA 

C,,  estimated 

a 

C,,  estimated 

11-20 

CMAE 

21-30 

CP 

C^  estimated 

P 

31-40 

GLD 

C  estimated 

h 

CARD  NO. 

11  (6F10.0) 

1-10 

CA 

linear  theory 

11-20 

CB 

roll  constants 

21-30 

CC 

31-40 

CD 

41-50 

XXREF 

51-60 

CNOT 

let  X  =  (X^  - 

XXREF) 

PHI  =  CA(X) 

+  CB(X^)  +  CC(X^) 

+  CD(X^)  +  CNOT 
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CARD  NO.  12  (7F10.0) 


1-10 

B(l) 

•^1 

11-20 

B(2) 

•^2 

21-30 

B(3) 

^1 

31-40 

B(4) 

^2 

41-50 

B(5) 

51-60 

B(6) 

"2 

61-70 

B(7) 

CARD  NO.  13 

(5F10.0) 

1-10 

B(8) 

“"2 

11-20 

B(9) 

21-30 

B(10) 

•^2 

31-40 

B(ll) 

41-50 

B(12) 

CARD  NO.  14 

(1311) 

NCON(J),  J  =  1,12, NEWS! 

CCON(J)  =  NCON(J)  J  =  1,12 

corresponds 

to  card  numbers  12  and  13 

NEWST  =  0 

program  will  restart  with 

NCON(ll)  and  NCON(12)  equal 

to  1 

NEWST  =  1 

no  restarts  allowed 

CARD  NO.  15 

(3011) 

NNCON(J)  J  =  1,30  for  ROLL-llKEVE 


ROLL-HELVE  INPUT 


CARD  NO.  1  (10A3) 
IBM  Col 


1-30 

SHOTN 

Title 

CARD  NO. 

2  (8F10.0) 

1-10 

DPT 

Time  Step 

11-20 

RHO 

Density 

21-30 

AS 

Speed  of  Sound 

31-40 

AMREF 

Reference  Mach  No, 

CARD  NO.  3  (8F1Q.0) 


1-10 

D1 

Reference  diameter,  inches 

11-20 

AIX 

2 

Axial  Moment  of  Inertia,  lb-in 

21-30 

AIY 

2 

Y-Transverse  Moment  of  Inertia,  lb-in 

31-40 

AIZ 

2 

Z-Transverse  Moment  of  Inertia,  lb-in 

41-50 

AIXY 

Product  of  Inertia 

51-60 

WGT 

Weight,  lbs. 

61-70 

ACG 

Offset  CG 

CARD  NO. 

4  (8F10.0) 

1-10 

XCGR 

Reference  CG,  percent  of  length 

11-20 

XLR 

Reference  length,  inches 

21-30 

XCG 

CG,  percent  of  length 

31-40 

XFIN 

No  of  Fins  or  Assymetries 

1~  .1^2^' 
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ROLL  HELVE  INPUT  (continued) 


C/'RD  NO.  5  (8F10.0) 


1-10 

TCON(l) 

,  degrees 

Q 

11-20 

TCON(2) 

q^,  deg/sec 

21-30 

ICON (3) 

,  degrees 

31-40 

ICON (4) 

LI 

r^,  deg/sec 

41-50 

TCON(5) 

C  - 
ma 

51-60 

TC0N(6) 

C  - 
mq 

61-70 

ICON (7) 

C  - 
npa 

71-80 

TCON(8) 

Not  used 

CARD  NO.  6^^  to 

6(Nt)  (8F10.0) 

Experimental  Data 

1-10 

TIME(K) 

Time,  sec 

11-20 

DIST(1,K) 

X  Distance,  ft 

21-30 

DIST(2,K) 

Y  Distance,  ft 

31-40 

DIST(3,K) 

Z  Distance,  ft 

41-50 

PHIX(K) 

0,  radians 

51-60 

THETFP(K) 

9pp,  degrees 

61-70 

PSIFP(K) 

Ypp,  degrees 

CARD  NO. 

1-10 

200. 

Stop  Card 

CARD  NO.  7,,..  to 

7(nnN) 

From  ROLL-NUT 

1-10 

XTIME(K) 

Time,  seconds 

11-20 

XQ(K) 

t;,  deg/ft 

7y 


roll  reeve  input  (continued) 


CARD  NO. 

to  7  ,  (8E10.0) 

^  (NNN)  _ 

(continued) 

21-30 

XR(K) 

R,  deg/ft 

31-40 

THE 

®RB’ 

41-50 

PRB 

'*'rb’ 

51-60 

XPHl(K) 

(8  radians 

G’ 

CARD  NO. 

1-10 

200. 

Stop  Gard 

CARD  NO.  8 

(8F10.0) 

1-10 

ACON(l) 

X  ,  feet 
o 

11-20 

AGON (2) 

X  ,  ft/sec 
o 

21-30 

AGON (3) 

Y  ,  feet 

0 

31-40 

AGON (4) 

Y  ,  ft/sec 
0 

41-50 

AGON (5) 

Z  ,  feet 
o 

51-60 

AGON (6) 

Z  ,  ft/sec 
o 

61-70 

AGON  (7).. 

s 

o 

71-80 

AGON (8) 

^Xa^ 

CARD  NO. 

q  (8F10.0) 

1-10 

AGON (9) 

^X62 

11-20 

AGON (10) 

in 

21-30 

AGON (11) 

S 

o 

31-40 

AC0N(12) 

^YP 

41-50 

AGON (13) 

51-60 

AGON (14) 

S 

o 

61-70 

AGON (15) 

^'Nn 

71-80 

AGON (16) 
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ROLL  HEEVE  INPUT  (conclud ed ) 


CARD  NO. 

10 

(8F10.0) 

1-10 

ACON(17) 

Cv 

Ypa 

11-20 

AC0N(18) 

SpB 

21-30 

ACON(19) 

31-40 

ACON(20) 

^Na 

41-50 

ACON(21) 

^Na^ 

51-60 

ACON(22) 

Spa 

61-70 

ACON(23) 

Not  Used 

71-80 

ACON(24) 

P  ,  rad/sec 
o 

CARD  NO. 

11 

(8F10.0) 

1-10 

ACON(25) 

0  ,  radians 

0 

11-20 

ACON(26) 

Sp 

21-30 

ACON(27) 

31-40 

ACON(28) 

A1-50 

AC0N(29) 

SPNFORM  (on 

51-60 

ACON{30) 

Not  Used 

CARD  NO. 

12 

(3011) 

1-30 

NNCON(l)  - 

NNCON(30) 
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R0LL-ANC.L1-:  INPUT 


CARD  NO.  1 

(10A3) 

1-30 

SHOTN 

CARD  NO.  2 

(1013) 

1-3 

LN 

A-6 

NA 

7-9 

ND 

10-12 

NB 

CARD  NO.  3 

(8F10.0) 

1-10 

DPT 

11-20 

RUO 

21-30 

AS 

31-AO 

AMREF 

CARD  NO.  4 

(8F10.0) 

1-10 

D1 

11-20 

AIX 

21-30 

AIY 

31-40 

AIZ 

41-50 

AIXY 

51-60 

WGT 

61-70 

ACG 

CARD  NO.  5 

(8F10.0) 

1-10 

XCGR 

11-20 

XLR 

21-30 

XCG 

31-40 

XFIN 

Title 

Starting  Pt.  of  fit 
Und  Pt.  of  Sectional  fit 
Incremental  section  interval 
Summing  section  interval 

Time  step 
De.asity 

Speed  of  sound 
Reference  Mach  No, 

Reference  diameter,  inches 
Axial  Moment  of  Inertia 
Y~Transverse  Moment  of  Inertia 
Z~Transverse  Momt’nt  of  Inertia 
Product  of  Inertia 
Weight,  lbs. 

Offset  CG 

Reference  CG,  percent  of  length 
Reference  lengtli,  indies 
CG,  percent  of  length 
No  of  Fins  or  Assymetries 
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ROLL  ANGLE  IIU’UT  (continued] 


CARD  NO.  6  (8F10.0) 

1-10 

TCON(l) 

,  degrees 

o 

11-20 

ICON (2) 

q^,  deg/sec 

21-30 

ICON (3) 

,  degrees 

o 

31-40 

ICON (4) 

r^,  deg/sec 

41-50 

ICON (5) 

C  - 
ma 

51-60 

ICON (6) 

c  - 

mq 

61-70 

ICON (7) 

c  - 

npa 

71-80 

ICON (8) 

Not  used 

CARD  NO.  7^^^  to  7^^^^  (8F10.0) 

1-10 

TIME(K) 

Time,  sec 

11-20 

DIST(1,K) 

X  Distance,  ft 

21-30 

DIST(2,K) 

Y  Distance,  ft 

31-40 

DIST(3,K) 

Z  Distance,  ft 

41-50 

PHIX(K) 

0,  radians 

51-60 

TllETFP(K) 

0pp.  degrees 

61-70 

PSIFP(K) 

Ypp,  degrees 

CARD  NO.  7^^,^^^^ 

1-10 

200. 

Stop  Card 

CARD  NO.  8  (ai’lO.O) 

1-10 

lI(:ON(i ) 

X  ,  feet 
o 

11-20 

HCON(2) 

X  ,  feet/sec 
o 

21-30 

HCON(3) 

Y  ,  feet 

0 

31-40 

HCON(4) 

Y  ,  feet/scc 
o 

ROLL  ANGLE  INPUT  (continued) 


CARD  NO.  8  (contd.)  (8F10.0) 


41-50 

HCON(5) 

t.  ,  feet 
o 

51-60 

HCON(6) 

Z^,  feet/sec 

61-70 

HCON(7) 

s 

o 

71-80 

HCON(8) 

CARD  NO. 

9  (8F10.0) 

1-10 

HCON(9) 

11-20 

HCON(IO) 

m 

21-30 

HCON(ll) 

0 

31-40 

HCON(12) 

Sa 

41-50 

HCON(13) 

51-60 

HCON(14) 

O 

61-70 

HCON(15) 

^Na 

71-80 

HCON(16) 

CARD  NO. 

10  (8fl0.0) 

1-10 

HCON(17) 

Spa 

11-20 

HCON(18) 

SpR 

21-30 

HCON(19) 

^Xa 

2 

31-40 

HCON(20) 

Sa 

41-50 

HCON(21) 

S 

51-60 

HCON(22) 

Cv  - 
ipa 

61-70 

11C0N(23) 

Not  Used 

71-80 

HCON(24) 

p  ,  rad/sec 

ROLL  ANGLE  INPUT  (continued) 


CARD  NO. 

11  (8F10.0) 

1-10 

HCON(25) 

0^,  radians 

11-20 

HCON(26) 

"IP 

21-30 

HCON(27) 

^16 

31-40 

HCON(28) 

^10 

41-50 

HCON(29) 

SPNFORM  (on  - 

51-60 

HCON(30) 

Not  used 

CARD  NO. 

12  (3011) 

1-30 

NHCON(l)  - 

NHCON(30) 

CARD  NO. 

“(1)  “  “(NT) 

1-10 

TIME(K) 

Time,  sec 

11-20 

VEL(K) 

Total  Velocity, 

21-30 

AQBAR(K) 

Dynamic  Pressure 

31-40 

TRB(K) 

®RB’  '^®8’'Res 

41-50 

PRB(K) 

degrees 

CARD  NO. 

13 

(NT+1) 

1-10 

200. 

Stop  Card 

CARD  NO. 

14  (8F10.0) 

1-10 

ACON(l) 

f‘RB  »  degrees 

11-20 

ACON(2) 

0 

deg/sec 

21-30 

ACON(3) 

^B  •  degrees 

31-40 

ACON(4) 

0 

r^,  deg/sec 
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ROLL  ANGLE  INPUT  (continued) 


CARD  NO.  lA  (cont.)  (8F10.0) 


41-50  ACON(5)  C 

mo 


51-60 

ACON(6) 

C 

ma 

61-70 

ACON(7) 

C 

ma^ 

71-89 

ACON(8) 

"no 

CARD  NO.  15 

(8F10.0J 

1-10 

ACON(9) 

"n3 

11-20 

ACON(IO) 

"n33 

21-30 

ACON(ll) 

C 

mq 

31-40 

ACON(12) 

41-50 

ACON(13) 

C 

nr 

51-60 

ACON(14) 

C 

nr2 

61-70 

ACON(15) 

C 

npa 

71-80 

ACON(16) 

C  „ 
mpB 

CARD  NO.  16 

(8F10.0) 

1-10 

ACON(17) 

"n0a 

11-20 

ACON(18) 

AIXY 

21-30 

ACON(19) 

AIX 

31-40 

ACON(20) 

41-50 

ACON(21) 

C 

mu  2 

51-60 

ACON(22) 

C 

na^ 

61-70 

AGON (23) 

"mo0 

71-80 

ACON(24) 
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ROLL  ANGLE  INPUT  (concluded) 


CARD  NO.  17 

(8F10.0) 

1-10 

ACON(25) 

Not  used 

11-20 

ACON(26) 

c  - 

ma 

21-30 

ACON(27) 

C  - 
ma 

3 

31-40 

ACON(28) 

C  - 
mq 

41-50 

ACON(29) 

C  - 
mq^ 

51-60 

ACON(30) 

c  - 

npa 

61-70 

ACON(31) 

Not  used 

71-80 

ACON(32) 

Not  used 

CARD  NO.  18 

(8F10.0) 

1-10 

ACON(33) 

Not  used 

11-20 

ACON(34) 

Not  used 

21-30 

ACON(35) 

Not  used 

31-40 

ACON(36) 

0  ,  radians 
o 

41-50 

ACON(37) 

n  ,  rad/sec 
o 

51-60 

ACON(38) 

^IP 

61-70 

ACON(39) 

^16 

71-80 

ACON(40) 

^10 

CARD  NO.  19 

(3011) 

1-40 

NNCON(l) 

-  NNCON(40) 

MANGLE  INPUT 


CARD  NO.  1  (II) 


1 

LMN 

No.  of  Rounds 

CARD  NO. 

2a 

(10A3) 

1-30 

SHOTN 

Title 

CARD  NO. 

3a 

(8F10.0) 

1-10 

DPT 

Time  step 

11-20 

RHO 

Density 

21-30 

AS 

Speed  of  sound 

31-40 

AMREF 

Reference  Mach  No. 

CARD  NO. 

4a 

(8F10.0) 

1-10 

D1 

Reference  diameter,  inches 

11-20 

AIX 

Axial  Moment  of  Inertia 

21-30 

AIY 

Y-Transverse  Moment  of  Inertia 

31-40 

AIZ  = 

Z-Transversc  Moment  of  Inertia 

41-50 

AIXY 

Product  of  Inertia 

51-60 

WGT 

Weight,  lbs 

61-70 

ACC 

Offset  CG 

CARD  NO. 

5a 

(8F10.0) 

1-10 

XCGR 

Reference  CG,  percent  of  length 

11-20 

XLR 

Reference  length,  inches 

21-30 

XCG 

CG,  percent  of  length 

31-40 

XFIN 

No.  of  Fins  or  Assymetries 
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MANGLE  INl’UT  (continued) 


CARD  NO.  6a  (8F10.0) 


1-10 

TCON(l) 

,  degrees 

o 

11-20 

ICON (2) 

q^,  deg/sec 

21-30 

TCON(3) 

,  degrees 

o 

31-40 

ICON (4) 

r  ,  deg/sec 
o 

41-50 

ICON (5) 

C  - 
ma 

51-60 

TCON(6) 

c  - 

mq 

61-70 

ICON (7) 

c  - 

npa 

71-80 

TCON(8) 

Not  used 

CARD  NO. 

7a^^  to  (8F10.0) 

1-10 

TIME(K) 

Time,  sec 

11-20 

DIST(1,K) 

X  Distance,  ft 

21-30 

DIST(2,K) 

Y  Distance,  ft 

31-40 

DIST(3,K) 

Z  Distance,  ft 

41-50 

PHIX(K) 

0,  radians 

51-60 

THETFP(K) 

Opp,  degrees 

61-70 

PSIFP(K) 

Ypp,  degrees 

CARD  NO. 

^^(NT+1) 

1-10 

200. 

Stop  Card 

CARD  NO. 

8a  (8F10.0) 

1-10 

HCON(l) 

X  ,  feet 
o 

11-20 

HCON(2) 

X^,  feet/sec 

21-30 

MCON(3) 

Y  ,  feet 
o 

Y  ,  feet/sec 
o 


31-40 


HCON(4) 


MANGLK  INl'UT  (continued) 


CARD  NO.  8a(contd.) 

(8F10.0) 

41-50 

HC0N(5) 

Z  ,  feet 
o 

51-60 

HCON(6) 

Z  ,  feet/ 
o 

61-70 

HC0N(7) 

'*0 

71-80 

I1C0N(8) 

^XOL^ 

CARD  NO.  9a  (8F10.0) 

1-10 

HCON(9) 

11-20 

IICON(IO) 

\ 

21-30 

IlCON(ll) 

31-40 

1IC0N(12) 

Sb 

41-50 

HCON(13) 

''YS3 

51-60 

11C0N(14) 

s 

61-70 

11C0N(15) 

0 

Sa 

71-80 

HCON(16) 

CARD  NO.  10  a  (8F10.0.) 

1-10 

I1C0N(17) 

11-20 

HCON(18) 

1  pn 

SpS 

21-30 

HCON(19) 

^Xa 

31-40 

HCON(20) 

2 

^Na 

41-50 

1IC0N(21) 

^Na 

51-60 

11C0N(22) 

S'pa 

61-70 

1IC0N(23) 

Not  U.sc'd 

71-80 

HCON(24) 

p^,  rad/se 
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MANGLE  INPUT  (continued) 


CARD  NO.  Ha  (8F10.0) 


1-10 

HCON(25) 

0^,  radians 

11-20 

HC0N(26) 

^IP 

21-30 

HCON(27) 

^16 

31-AO 

HCON(28) 

^10  I  _  I 

Al-50 

11C0N(29) 

SrNFORM  (on  - -) 

X 

51-60 

IICON(JO) 

Not  used 

CARD  NO.  12a  (3011) 

1-30 

NllCON(l)  -  NHCON(30) 

CARD  NO.  13a^^^  to  13a, (8F10.0) 

1*-10 

TIME(K) 

Time,  sec 

11-20 

VEL(K) 

Total  Velocity,  ft/sec 

21-30 

AQBAR(K) 

Dynamic  Pressure 

31-AO 

TRBtK) 

®RB’  '^^Si'ees 

Al-50 

PRH(K) 

degrees 

CIVRD 

1-10  200.  Stop  Card 


CARD  NO.  14a  (8F10.0) 


1-10 

ACON(l) 

®Rli  ’ 
o 

o 

1 

A(:ON(2) 

,,  ,  deg/sec 

0 

21-30 

ACON(3) 

,  degrees 

31-AO 

ACON(4) 

VI 

r  ,  deg/sec 

O 

MANGLE  IKl’UT  (continued) 


CARD  NO. 

14a(cont.)  (SEIO.O) 

41-50 

ACON(5) 

C 

mo 

51-60 

AGON (6) 

C 

ma 

61-70 

ACON(7) 

C 

"’«3 

71-80 

AGON (8) 

^no 

CARD  NO. 

15a(8F10.0) 

1-10 

AGON (9) 

11-20 

ACON(IO) 

S03 

21-30 

ACON(ll) 

C 

mq 

31-40 

ACON(12) 

41-50 

ACON(13) 

C 

nr 

51-60 

AGON (14) 

C 

""2 

61-70 

AGON (15) 

C 

npa 

71-80 

ACON(16) 

C  „ 
mp6 

CARD  NO. 

16a(8F10.0) 

1-10 

ACON(l/) 

11-20 

ACON(18) 

AIXY 

21-30 

AC0N(19) 

AIX 

31-40 

ACON(20) 

41-50 

ACON(21) 

C 

ma^ 

51-60 

ACON(22) 

C 

na^ 

61-70 

ACON(23) 

^mot) 

71-80 

AGON (24) 

^mo0 
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MANGLE  INPUT  (concluded) 


CARD  NO.  17a 

(8F10. 

91 

1-10 

AC0N(25) 

Not  used 

11-20 

ACON(26) 

21-30 

ACON(27) 

c  - 

mot^ 

31-AO 

ACON(28) 

C  - 
mq 

41-50 

ACON(29) 

C  - 

mq2 

51-60 

ACON(30) 

c  - 

npa 

61-70 

ACON(31) 

Not  used 

71-80 

ACON(32) 

Not  used 

CARD  NO.  18a 

(8F10. 

01 

1-10 

AGON (33) 

Not  used 

11-20 

ACON(34) 

Not  used 

21-30 

ACON(35) 

Not  used 

31-40 

ACON(36) 

0  ,  radians 
o 

41-50 

ACON(37) 

p^,  rad/sec 

51-60 

ACON(38) 

Oi 

Ip 

61-70 

ACON(39) 

^16 

71-80 

ACON(40) 

^10 

CARD  NO.  19a 

(4011) 

1-40 

NNCON(LM,l) 

-  NNCON(LM,40) 
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Cards  2  through  19  must  be  read  'UIN' 
(i.e.  2a-19a,  2b-19b,  ....)• 


times,  stacked  behind  each  other 


cards  8-12  are  read  to  get  the  IICONS  to  bo  used  by  MANCIJi. 

cards  14-19  are  nert  to  supply  the  co,™.o„  coelflclents  for  the  multiple 
fit  aud  also  the  unlquo  coefficients  (ACON's  1.  2.  3,  4,  5.  8.  18,  19,  36, 

37)  of  the  first  data  set. 

►  A  fr,r  remaining  data  sets  to  supply  the 

cards  20-22  (Below)  are  repeated  for  the  remaining 

unique  coefficients  for  those  rounds. 


CARD  NO. 

20  (sno.o) 

1-10 

ACON(l) 

®RB  ’ 

11-20 

AC0N(2) 

q  ,  deg/sec 
o 

21-30 

AC0N(3) 

VI'  ,  degrees 

^  RB  ’ 
o 

31-AO 

ACON(A) 

r  ,  deg/sec 

0 

41-50 

AC0N(5) 

C 

mo 

51-60 

AC0N(8) 

C 

no 

61-70 

AC0N(18) 

AIXY 

71-80 

AC0N(19) 

AIX 

CARD  NO.  21  (8F10.01 


1-10 


AC0N(36) 


0  ,  radians 
o 

n  ,  rad/sec 
'  o 


11-20 


AGON (37) 


EXA>a’LE  DATA  DECK  Sl-TUr 

TWO  ROUNDS 

Card  1  (LMN  =  2) 

Cards  2a-19a 

Written  by  ROLL-HEAVE 

Cards  2b-19b 

Cards  8-12  HCONS  for  MANGLE 

Cards  14-19 

ACONS  for  MANGLE 

Cards  20-22 

THREE  ROUNDS 

Card  1  (LMN  =  3) 

Cards  2a-19a 

Cards  2b-19b  Written  by  ROLL-HEAVE 

Cards  2c-19c 

Cards  8-12  HCONS  for  MANGLE 

Cards  14-19 

Cards  20-22  ACONS  for  MANGLE 

Cards  20-22 
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MUL-HEEV  INPUT 


IBM  Col. 

CARD  NO. 

1  (11) 

1 

LMN 

Number  of  Rounds 

CARD  NO. 

_2a 

1-30 

SHOTN 

Title 

CARD  NO. 

3a 

(8F10.0) 

1-10 

DPT 

Time  Step 

11-20 

RHO 

Density 

21-30 

AS 

Speed  of  Sound 

31-40 

AMREF 

Reference  Mach  No. 

CARD  NO. 

4a 

(8F10.0) 

1-10 

D1 

Reference  diameter,  inches 

11-20 

AIX 

2 

Axial  Moment  of  Inertia,  Ib-ln 

21-30 

AlY 

Y-Transverse  Moment  of  Inertia, 

31-40 

AIZ 

Z-Transverse  Moment  of  Inertia, 

41-50 

AIXY 

Product  of  Inertia 

51-60 

WGT 

Weight,  lbs 

61-70 

ACG 

Offset,  CG 

CARD  NO. 

5a 

(8F10.0) 

1-10 

XGGR 

Reference  CG,  percent  of  length 

11-20 

XLR 

Reference  length,  Indies 

21-30 

XCG 

CG,  percent  of  length 

31-40 

XFIN 

No.  of  Fins  or  Assymetries 
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CARD  NO. 6a 

(8K10.0) 

1-10 

TCON(l) 

®R3  >  degrees 

11-20 

ICON (2) 

o 

Rq.  deg/sec 

21-30 

ICON (3) 

'*'33  .  degrees 

31-AO 

ICON (4) 

0 

deg/sec 

41-50 

ICON ( 5 ) 

C  - 

ma 

51-60 

ICON (6) 

C  - 

mq 

61-70 

3CON(7) 

C  - 

npa 

71-80 

're  UN  (8) 

Not  used 

CARD  NO.  7;i 


to  7a 


"'(]n)  Experimental  Data 


1-10 

TIME(K) 

Time,  sec 

11-20 

DIST(1,K) 

X  Distance,  ft 

21-30 

D1ST(2,K) 

Y  Distance,  ft 

31-40 

D1ST(3,K) 

Z  Distance,  ft 

41-50 

PHIX(K) 

0,  radians 

51-60 

Tl(ETFP(K) 

Oj,.p,  degrees 

61-70 

P.Siri’(K) 

degrees 

CARD  NO. 

7‘D 

(NT+1) 

1-10 

200. 

5Lo|)  Card 

CARD  NO. 

Prom  KOLi.-NUT 

1-10 

XTIME(K) 

Time,  sec 

11-20 

XQ(K) 

Q,  deg/ft 

CARD  NO.  8n  to  8a 

( K  )  1 

^fUL-lll•EV 

INl'UT  (continued) 

(continued) 

21-30 

XR(K) 

R,  deg/ft 

31-40 

TR» 

degrees 

41-50 

rRi» 

'•'Rjjt  degrees 

51-60 

XPIII(K) 

0^,  radians 

CARD  NO.  8a, 

(NNN+l) 

1-10 

200. 

Stop  Card 

CARD  NO.  9a  (8F10.0) 

1-10 

ACON(l) 

feet 

11-20 

ACON(2) 

X^,  ft/sec 

21-30 

ACON(3) 

Y  ,  feet 
o 

31-40 

A(:ON(4) 

V^,  ft/sec 

41-50 

ACON(5) 

2^,  feet 

51-60 

AC0N(6) 

ft/sec 

61-70 

ACON(7) 

71-80 

AC0N(8) 

O 

CARD  NO.  lOu  (8F10.0) 

1-10 

ACON(9) 

^•X8 

11-20 

ACON(IO) 

L 

21-30 

AGON (11) 

ni 

s- 

31-40 

AGON (12) 

0 

U-50 

AGON (13) 

51-60 

agon (14) 

3 

^'N 

51-70 

AGON (15) 

0 

G 

No 

71-80 

AGON (16) 

G 

MUL-HEEV  INrUT  (continued) 


CARD  NO. 

11a 

(8F10.0) 

1-10 

AC0N(17) 

Spa 

11-20 

AC0N(18) 

SpS 

21-30 

AC0N(19) 

Sa2 

31-40 

ACON(20) 

^Na 

41-50 

AC0N(21) 

^Na^ 

51-60 

ACON(22) 

Spo 

61-70 

ACON(23) 

Not  Used 

71-80 

ACON(24) 

P  ,  rad/sec 
o 

CARD  NO. 

12a 

(8F10.0) 

1-10 

ACON(25) 

0  ,  radians 
o 

11-20 

ACON(26) 

^IP 

21-30 

ACON(27) 

31-40 

ACON(28) 

^10 

41-50 

ACON(29) 

SPNFORM  (on 

51-60 

AC0N(30) 

Not  used 

CARD  NO. 

13a 

(3011) 

1-30 

NNC0N(LM,1) 

-  NNCON(LM,30) 

o  Cards  2  through  13  must  be  read  'LMN'  times  stacked  behind  each  other 
(i.e.,  2a-13a,  2b-13b,  ...). 

o  Cards  9  through  13  are  then  read  once  again  containing  the  initial  estimates 
for  the  common  aerodynamics  and  for  the  unique  coefficients  (ACON(l),  (2), 
(3),  (4),  (5),  (6),  (11),  (14),  (24),  (25),  (29))  pertaining  to  tlic  first 
data  set. 
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Cards  14a,  15a,  16a  are  read  for  the  unique  coefficients  and  NNCON's  for 
the  remaining  data  sets. 


CARD  NO.  14a 

(8F10.0) 

1-10 

ACON(l) 

X  ,  feet 

0 

11-20 

ACON(2) 

X  ,  ft/sec 
o 

21-30 

AC0N(3) 

Y  ,  feet 

0 

31-40 

ACON(4) 

Y  ,  ft/sec 
0 

41-50 

ACON(5) 

Z  ,  feet 
o 

51-60 

ACON(6) 

Z  ,  ft/sec 

0 

61-70 

ACON(ll) 

S 

0 

71-80 

AC0N(14) 

s 

0 

CARD  NO.  15a 

(8F10.0) 

1-10 

ACON(24) 

P  ,  rad/sec 
o 

11-20 

ACON(25) 

0^,  radians 

21-30 

ACON(29) 

SPNFORll  (on 

CARD  NO.  16a 

(1111) 

1-11 

NNC0N(LM,1) 

-  NNCON(LM,30) 

I 


MUL-HEEV  INPUT  (Concluded) 


EXAMPLE  or  DATA  DECK  SETUPS 


TWO  ROUNDS 

Card  No.  1  LMN  =  2 

Card  No.  2a-13a 
Card  No.  2b-13b 

Card  No.  9a-13a 
Card  No.  14a-16a 

THREE  ROUNDS 

Card  No.  1  LMN  =  3 

Card  No.  2a-13a 
Card  No.  2b-13b 
Card  No.  2c-13c 

Card  No.  9a-13a 
Card  No.  14a-16a 


Written  by  ROLANCLE 


Coefficient  for  MUL-HEEV 


Written  by  ROLANCLE 


Card  No.  14b-16b 


LIST  01  SYMBOLS 


*0’  V" 

hnf'  -'mr'  *^iiir 


^n'  'nr  '^m 


dt 

in 


V 


CO  "p 

r 


^'RB'  '^RB’  '^RB 
p,  q,  r 


xy 


I’ 


L 


-  unit  vectors  aligned  with  the  earth  coordinate  system 

-  unit  vectors  aligned  with  the  fixed  plane  coordinate  system 
unit  vectors  aligned  with  tlie  body  fixed  coordinate  system 

-  derivative  witli  respect  to  time  in  body  lixed  coordinates 

-  mass  of  the  missile 

-  total  velocity  vector  relative  to  earth 
total  angular  velocity  relative  to  earth 

-  total  angular  momentum  of  the  missile 

-  body  fixed  velocity  components  relative  to  earth 

-  body  fixed  angular  velocity  components  relative  to  earth 

-  mass  moment  of  inertia  about  the  x-axis 

-  mass  moment  of  inertia  about  the  y-axis 

-  mass  moment  of  inertia  about  the  z-axis 

-  cross  product  of  mass  moment  of  inertia 

-  indicates  the  derivative  of  u  with  respect  to  time 

-  body  fixed  total  lorce  vector 
body  fixed  total  moment  vector 
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Fx,  Fy,  F, 

Lp,  Lq,  L^ 

®FP’  "^FP 

0 

'^‘m’  '^nr  "'m 


q 

A 

d 

0' 

t 

P 


^RB’  '^RB 

V  V  V 

Vx>  ’'yi 


—  body  fixed  force  components 

—  body  fixed  moment  components 

—  fixed  plane  Euler  angles 

—  roll  angle  about  the  body  axis 

—  missile  fixed  velocity  components  relative  to  the  trajectory 

—  dynamic  pressure 

—  reference  area  of  the  body 

—  reference  diameter  of  the  body 

—  aerodynamic  roll  angle 

—  time 

—  air  density 

—  body  fixed  Euler  angles 

—  Earth  fixed  velocity  components 


: 

: 


f  j 

P 

i 

I 
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INITIAL  DISTRIBUTION 


Hq  USAF/RDQRNl 

Hq  USAF/SAMI 

Hq  USAF/XOXFCM 

AFSC/IGFG 

AFSC/SDZA 

Hq  AFSC/DLCAW 

AFML/DO/ANIIC 

Hq  4950  TW/TZHM 

AFIT/LD 

ASD/ENYS 

,\SD/ENAZ 

AFFDL/FES 

TAC/DRA 

SAC/LGWC 

Hq  SAC/NRI ,  STINFO  Lib 
WRAMA/NIMEBL 
CRE/ADD/Publ ications 
AFWL/LR 

AUL  (ALU/LSE-70-239) 

SAPRI-LW-A 

AMXS\-DD 

AMXSY-A 

DRXBR-TE 

Lib,  K2400 

SARPA-TS-TS 

Dahlgren  Laboratory 

White  Oak  Laboratory 

NAV  ORD  STN/Tech  Lib 

NAV  WEAPONS  STN/20323 

NAV  Underwater  Systems  Ctr/Code  54 

USN  WEA  CNTR/Code  233 

Naval  Weapons  Center/Code  31 

Air  Force  Wpns  Lab/Tech  Lib 

NAV  AIR  SYS  COMD/Code  AIR- 532 

Office  Naval  Research/Code  473 

NASA  STINFO  FAC 

Lawrence  Rad  Lab/Dept  L-3 

The  John  Hopkins  Univ/ 

Applied  Physics  Lab 
Battelle  Memorial  Inst/Reports  Lib 
Inst  for  Defense  Analysis/ 
Classified  Lib 
Sandia  Laboratories 
The  RAND  Corp/Lib-D 
DDC/TC 
USAFTFWC/TA 
SARWV-RDT-L 


1  Naval  Weapons  Center/Code  3123  1 

1  Ogden  ALC/Nf'IWM  2 

1  AF  Spec  Comm  Cntr/SUR  3 

1  DAMA-WSA  1 

1  SARPA-FR-S-A  1 

1  US  Atomic  Energy  Comm/Hq  Lib  1 

1  AEDC/ARO,  Inc/DLCS  1 

1  AMXSY/DS  1 

1  AMCRD/WM  1 

1  Nav  t^ns  Eval  Fac/Code  WT  1 

1  Office  of  the  Chief  of  Nav  Opns/  1 

1  (OP-982E) 

1  Naval  Research  Lab/ Code  2627  1 

1  Hq  PACAF/LGWSE  4 

1  USAFTAWC/TE  1 

1  TAWC/TRADOCLO  1 

2  AFATL/DL  1 

1  AFATL/DLY  1 

1  AFATL/DLOU  1 

1  AFATL/DLOSL  9 

1  AFATL/DLYV  1 

1  AFATL/DLDL  1 

1  AFATL/DLDA  1 

1  AFATL/DLDE  1 

1  AFATL/DLDT  1 

1  AFATL/DLDG  1 

2  AFIS/INTA  1 

1  Naval  Sea  Systems  Comd/Code  SEA-  1 

1  0332) 

1  Naval  Sea  Systems  Comd/Code  SEA-  1 

2  992E) 

1  Nav  Ord  Lab/White  Oak  2 

1  Naval  Wpns  Center/Code  32602  1 

1  Naval  W^ns  Center/Code  3163  1 

1  Ogden  ALC/NMWRA  2 

1  AFLC/NMWMC  1 

1  ASD/ENESS  1 

1  AFATL/DU  2 

AFATL/DLA  1 

1  ADTC/SDC  1 

1  Redstone  Sci  Info  Cntr/Doc  Sec  1 


General  Electric  Co /Armament  Dept  10 

2 

1 

2 

1 

1 


105 


